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Chapitre  1:  Généralités  sur  les  modes  de  transfert  de  chaleur 


1.  Introduction: 

Parmi  les  modes  d’échange  d’énergie,  le  transfert  de  chaleur  est  certainement  l’un  des  plus  courants.  11 
intervient  naturellement  entre  deux  systèmes  dès  qu’une  différence  de  température  apparaît  entre  eux 
et  cela  quel  que  soit  le  milieu,  même  vide,  qui  les  sépare.  L’importance  de  ces  phénomènes  de 
transmission  de  chaleur  résulte  de  leur  intervention  dans  presque  tous  les  procédés  des  différents 
secteurs  d’activités:  chimie,  pharmacie,  agro-alimentaire,  environnement,  pétrole... Le  rôle  de 
transfert  de  chaleur  devient  souvent  essentiel  et  même  déterminant  lorsqu’il  est  à l’origine  des 
techniques  utilisées  comme,  par  exemple,  les  échangeurs  de  chaleur,  l’isolation  thermique  ou  le 
calorifugeage. . . 

On  distingue  trois  modes  de  transfert  de  chaleur:  conduction , convection  et  rayonnement. 

Dans  la  plupart  des  situations  rencontrées  dans  l’industrie,  ces  différents  modes  ne  sont  pas  dissociés 
et  peuvent  intervenir  ensemble. 

► Exemples: 

► Echange  de  chaleur  entre  deux  fluides,  séparés  par  une  paroi  solide,  circulant  dans  un 
échangeur  de  chaleur. 

► Transmission  de  la  chaleur  des  gaz  chauds  (produits  de  combustion)  vers  les  parois  internes 
d’une  chambre  de  combustion  d’une  chaudière  ou  d’un  four. 

► Déperdition  de  chaleur  dans  des  conduites  transportant  des  fluides  caloporteurs. 

► Déperdition  de  chaleur  dans  un  bâtiment  à travers  les  murs  en  contact  avec  l’extérieur. 

► Diffusion  de  chaleur  dans  un  local  via  un  radiateur  ou  un  plancher  chauffant. 

Lorsque  l’un  des  modes  est  prédominant  on  négligera  évidemment  les  deux  autres. 

Tous  ces  modes  de  transfert  de  chaleur  ont  la  même  conséquence:  la  variation  de  la  température  dans 
le  milieu  matériel  où  ils  interviennent.  Parfois,  cette  variation  de  température  peut  conduire  à un 
changement  de  phase  (vaporisation,  condensation,  fusion,  solidification).  Dans  cette  situation,  la 
chaleur  latente  joue  le  rôle  d’une  source  de  chaleur  (positive  ou  négative)  supplémentaire  qui  accroît  la 
complexité  du  phénomène  de  transfert  de  chaleur.  11  faut  noter  que  l’intensité  de  ce  transfert  devient 
beaucoup  plus  grande  en  présence  de  changement  de  phase,  ce  qui  explique  l’utilisation  de  ce 
phénomène  dans  de  nombreuses  applications  pratiques  utilisant  des  technologies  modernes  mettant  en 
jeu  des  flux  de  chaleurs  très  élevées. 

2.  Grandeurs  thermiques: 

2.1.  Chaleur: 

La  chaleur  est  une  forme  d’énergie  dont  les  manifestations  ont  la  particularité  d’être  directement 
ressenties  par  le  corps  humain.  Les  notions  de  chaud  ou  de  froid  résultent  de  la  perception  par  nos  sens 
de  l’accroissement  ou  la  diminution  de  la  chaleur. 

On  distingue  deux  formes  de  chaleur:  chaleur  sensible  et  chaleur  latente. 

2.1.1.  Chaleur  sensible: 

On  parle  de  chaleur  sensible  lorsqu’un  corps  qui  reçoit  ou  cède  de  la  chaleur  s’échauffe  ou  se  refroidit 
sans  changer  d’état. 

Cette  quantité  de  chaleur  est  donnée  par: 

Q = m x cp  x AT  (1) 

où  m est  la  masse  du  coips,  cp  sa  chaleur  massique  et  AT  la  variation  de  sa  température. 
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2.1.2.  Chaleur  latente: 

On  parle  de  chaleur  latente  lorsqu’un  coips  qui  reçoit  ou  cède  de  la  chaleur  l’utilise  pour  changer 
d’état  à température  constante. 

Cette  quantité  de  chaleur  est  donnée  par: 

Q = m x L (2) 


où  L est  le  coefficient  de  chaleur  latente. 

2.2.  Température: 

D’une  manière  générale,  la  température  est  une  mesure  de  l’état  de  vibration  (atomes)  ou  d’agitation 
(molécules)  des  constituants  fondamentaux  de  la  matière. 

En  effet,  l’énergie  thermique  d’un  milieu  matériel  correspond  à l’énergie  cinétique  de  ses  constituants 
fondamentaux  ayant  une  certaine  liberté  de  mouvement  (molécules,  atomes,  électrons  libres,  etc.). 
Ceux-ci  pourront  échanger  toute  ou  partie  de  leur  énergie  thermique,  c'est-à-dire  gagner  ou  perdre  de 
l’énergie  cinétique. 

La  température,  notée  T,  s’exprime  en  °C  ou  en  K. 

2.2.1.  Champ  de  température: 

Les  transferts  de  chaleur  sont  déterminés  à partir  de  l’évolution,  dans  l’espace  et  dans  le  temps,  de  la 
température.  Ainsi,  on  peut  définir  en  chaque  point  M d’un  coips  (solide,  liquide  ou  gaz)  un  scalaire, 
fonction  des  coordonnées  du  point  M et  du  temps,  appelé  champ  de  température:  T(M,  t). 

On  distingue  deux  cas  : 

► Champ  de  température  indépendant  du  temps  : le  régime  est  dit  permanent  ou  stationnaire  ; 

► Evolution  du  champ  de  température  avec  le  temps  : le  régime  est  dit  variable  ou 
instationnaire. 

2.2.2.  Surfaces  isothermes: 

On  appelle  surface  isotherme  le  lieu  des  points  ayant  à chaque  instant  la  même  température.  En  régime 
permanent,  les  isothermes  sont  invariantes',  en  régime  variable,  elles  sont  mobiles  et  déformables. 

Le  vecteur  gradient  de  température  est  en  tout  point  normal  à la  surface  isotherme  passant  par  ce  point. 


2.3.  Flux  de  chaleur: 

Pour  tous  les  modes  de  transfert  de  chaleur,  on  définit  le  flux  de  chaleur  ou  la  puissance  thermique, 
noté  0,  comme  la  quantité  de  chaleur  Q transmise  à travers  une  surface  d’aire  S pendant  l’intervalle 
de  temps  A t : 

(HO  (3) 

Comme  le  flux  de  chaleur  n’est  pas  nécessairement  uniforme  sur  toute  la  surface  d’échange  5,  on 
introduit  la  notion  de  densité  de  flux  local  (p  qui  est  défini  comme  la  puissance  thermique  traversant 
l’unité  de  surface: 
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<P=j  ( W/m 2)  (4) 

3.  Aperçu  des  modes  de  transfert  de  chaleur: 

3.1.  Conduction: 

On  sait  que  l’énergie  interne  d’un  coips  solide,  liquide  ou  gazeux  correspond  à l’énergie  cinétique  de 
ses  éléments  constituants  qui  se  trouvent  dans  un  mouvement  permanent  (atomes,  molécules, 
électrons).  La  conduction  représente  le  processus  de  transmission  de  chaleur  par  le  contact  direct  entre 
les  particules  d’un  corps  donné  (qu’il  soit  solide,  liquide  ou  gazeux)  ou  entre  des  coips  ayant  des 
niveaux  énergétiques  (donc  des  températures)  différents. 

► Exemples-. 

► Transmission  de  la  chaleur  à travers  une  paroi  d’échange  séparant  deux  fluides  circulant  dans 
un  échangeur  de  chaleur. 

► Transfert  de  chaleur  à travers  les  parois  d’un  bâtiment  ou  des  conduites  transportant  des 
fluides  chauds  ou  froids. 

Dans  n’importe  quel  milieu  matériel,  ce  mode  de  transfert  de  chaleur  tend  à uniformiser  la 
température,  donc  la  répartition  d’énergie  cinétique  des  particules  constituantes  dans  la  masse  du 
milieu. 

► Remarques'. 

► La  conduction  est  le  seul  mode  à exister  dans  un  solide  opaque-,  ainsi  la  conduction  concerne 
essentiellement  les  solides.  Exemple : transmission  de  la  chaleur  à travers  une  paroi  en  métal 
ou  un  mur  en  béton.  Elle  est  aussi  prédominante  dans  les  fluides  au  repos  (la  convection  est 
totalement  exclue). 

► Pour  les  solides  transparents , une  partie  de  la  chaleur  peut  être  transmise  par  rayonnement. 
Exemple:  transmission  de  la  chaleur  à travers  un  vitrage  (porte  ou  fenêtre  vitrée). 

► Dans  les  liquides  et  les  gaz , la  convection  et  le  rayonnement  peuvent  se  superposer  à la 
conduction  qui  est  négligeable.  Exemple : transmission  de  la  chaleur  à travers  les  gaz  chauds 
dans  une  chambre  de  combustion  ou  dans  un  four. 

► Flux  de  chaleur  conductif:  Loi  de  Fourier 

La  conduction  thermique  est  régie  par  une  loi  dite  de  Fourier  qui  exprime,  dans  le  cas  de  milieux 
homogènes  et  isotropes,  la  proportionnalité  du  vecteur  densité  surfacique  de  flux  de  chaleur  et  du 
gradient  local  de  température  grad(T ): 


ç{M)  = -Ax  grad(T(M))  ( W/m 2)  (5) 

Le  coefficient  de  proportionnalité  A est  une  caractéristique  physico-chimique  du  milieu  matériel 
désignée  sous  le  nom  conductivité  thermique.  Elle  s’exprime  en  W/m.°C  ou  en  Kcal/h.m.°C  et 
traduit  l’aptitude  des  matériaux  à conduire  la  chaleur. 

La  densité  de  flux  de  chaleur  <px  dans  une  direction  caractérisée  par  un  vecteur  unitaire  ex  est: 


çx(M)  = -Ax  grad(T(M)).ex  = -A  x — (M)  (6) 

Le  signe  (-)  est  lié  au  fait  que  le  flux  de  chaleur  est  positif  quand  le  gradient  thermique  est  négatif.  La 
chaleur  s’écoule  de  la  zone  la  plus  chaude  vers  la  zone  la  plus  froide. 
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La  conductivité  thermique  A dépend  principalement  de  la  nature  du  matériau.  D’une  manière  générale, 
les  métaux  sont  beaucoup  plus  conducteurs  de  la  chaleur  que  les  coips  non  métalliques.  Les  gaz  sont 
de  mauvais  conducteurs. 


La  conductivité  thermique  dépend  également  de  la  température  et  de  l’état  physique  du  matériau 
(solide,  liquide,  gaz).  Par  exemple,  la  conductivité  de  l’eau  à: 

► l’état  gazeux  (vapeur)  est:  A~  0,  025  W /m.  °C  (à  T = 100°C), 

► l’état  liquide : A ~ 0,  6 W/m.  °C  (à  T = 20 °C), 

* l’état  solide  (glace)  est:  A~2,2W/m.°C(à.T  = 0 °C). 

11  existe  d’autres  facteurs  qui  influencent  la  conductivité  comme  l’humidité  et  la  densité. 

Le  flux  de  chaleur,  noté  0,  représente  la  quantité  de  chaleur  qui  passe  par  unité  de  temps  à travers  une 
surface  isotherme  5. 11  s’exprime  par  la  relation  : 


0 = (p  x dS  = (p  x S = -S  x A x — ( W ) (7) 

Js  dx 

> Exercice  d’application  1: 

On  considère  un  four  industriel  dont  les  parois  sont  construites  en  briques  réfractaires  d’épaisseur 
e = 15  cm  et  de  conductivité  thermique  A = 1, 7 W/m.  °C.  Des  mesures  réalisées  lors  du 
fonctionnement  du  four,  en  régime  permanent,  ont  révélé  les  températures  suivantes:  1127°C et 
877°C  respectivement  sur  les  surface  frontières  interne  et  externe  d’une  paroi  du  four.  Calculer  le 
flux  de  chaleur  perdu  à travers  la  paroi  du  four  de  0,  5 m de  hauteur  et  1,  2 m de  longueur. 


Solution: 

• Paroi  simple  en  briques  réfractaires  (X  = 1,  7 W/m.  °C , e = 0, 15  m,  L = 0,  5 m et  H = 

1, 2 m) 

■ Surfaces  frontières  de  la  paroi  simple:  ( Si=Se  = S = LxH , T Pi  = 1127°C  et  T Pe  = 
877°C) 


La  paroi  en  briques  réfractaires  correspond  à un  milieu  matériel  conductif  sans  sources  ni  puits  de 
chaleur,  en  régime  permanent,  pour  lequel  la  chaleur  se  transmet  suivant  une  seule  direction 
( e « L et  H).  Le  flux  de  chaleur  traversant  cette  paroi  s’exprime  comme  suit  : 
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0=Â.S.— — — = A.H.L.— — — * 1700  W = 1,7  KW 
e e 

Les  pertes  thermiques  à travers  cette  paroi  sont  élevées,  on  devra  donc  penser  à réaliser  une  isolation 
thermique  efficace  de  cette  paroi  ( voir  chapitre  3). 

3.2.  Convection: 

La  convection  est  le  processus  de  transmission  de  chaleur  qui  s’appuie  sur  un  milieu  matériel,  avec 
mouvement  macroscopique  de  matière.  11  ne  peut  donc  avoir  lieu  qu’au  sein  des  liquides  ou  des  gaz; 
c’est  le  fluide  en  mouvement  qui  transporte  la  chaleur. 

La  convection  peut  apparaître  entre  des  couches  de  fluide  ayant  des  températures  différentes  ou  entre 
un  fluide  en  écoulement  au  long  d’une  paroi  solide  ayant  une  température  différente  de  celle  du  fluide. 
La  convection  est  le  mécanisme  le  plus  important  de  transfert  de  chaleur  entre  une  paroi  solide  et  un 
fluide. 

On  distingue  traditionnellement  trois  types  de  convection: 

► La  convection  forcée  qui  se  produit  quand  le  mouvement  de  fluide  est  provoqué  par  une 
action  externe  telle  qu’une  pompe  ou  un  ventilateur.  C’est  le  cas  des  échangeurs  de  chaleur 
industriels,  des  radiateurs  de  voiture  qui  sont  essentiellement  des  convecteurs:  deux  fluides  en 
circulation  forcée  échangent  de  la  chaleur  à travers  une  paroi  solide  dont  la  température  diffère 
des  leurs.  Comme  les  vitesses  d’écoulement  peuvent  atteindre  des  valeurs  très  élevées,  le 
transfert  par  convection  forcée  est  souven  t extrêmemen  t efficace. 

► La  convection  naturelle  ou  libre  qui  apparaît  spontanément  sous  l’influence  de  différences  de 
densités  dues  à des  différences  de  températures  au  sein  du  fluide  situé  dans  le  champ  de 
gravité.  Le  fluide  chaud,  plus  léger,  tend  alors  à montrer,  tandis  que  le  fluide  froid,  plus  dense, 
descend.  La  configuration  et  l’intensité  de  la  convection  naturelle  sont  en  relation  directe  avec 
les  conditions  thermiques  qui  la  déclenchent,  la  nature  du  fluide  et  le  volume  de  l’espace  où  a 
lieu  ce  processus.  Le  mécanisme  de  la  convection  naturelle  se  produit  dans  plusieurs 
situations.  11  est,  par  exemple,  responsable  de  l’homogénéisation  de  la  température  dans  un 
local  d’habitation:  transfert  de  chaleur  du  radiateur  vers  les  zones  froides  du  local. 

► La  convection  mixte  qui  correspond  au  couplage  des  deux  phénomènes  précédents 
(convections  naturelle  et  forcée)  quand  les  vitesses  d’écoulement  fictives,  dues  aux  deux  types 
de  convection  sont,  considérées  séparément,  du  même  ordre  de  grandeur. 

► Flux  de  chaleur  convectif:  Loi  de  Newton 

Le  flux  de  chaleur  échangé  par  convection  entre  une  paroi  solide  (de  surface  5)  à température  Ts  et  un 
fluide  à température  Tœ  (température  mesurée  loin  de  la  paroi)  peut  être  exprimé  au  moyen  de  la 
relation  suivante: 

0=hxSx(Ts  — Tœ)  (8)  ( Loi  de  Newton ) 

où  h est  le  coefficient  d’échange  de  chaleur  par  convection  (unité:  W/m2.°C  ou  K cal /h.  m2.  °C ). 

Le  coefficient  h dépend,  en  général,  de  l’état  (lisse  ou  rugueux)  et  de  la  forme  géométrique  de  la  paroi 
solide,  de  la  nature  de  l’écoulement  (convection  forcée  ou  naturelle,  laminaire  ou  turbulente)  et 
essentiellement  des  propriétés  thermophysiques  du  fluide  {p,  cp,  A,  ju)  et  de  l’écart  de  température 
(AT  = TS-Tœ). 


On  donne  à titre  indicatif,  dans  le  tableau  suivant  les  ordres  de  grandeur  des  valeurs  atteintes  par  le 
coefficient  h pour  différentes  configurations. 
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Tableau: 


Configuration 

h (W/m2.°C) 

Convection  naturelle 

Dans  les  gaz 

2-30 

Dans  les  liquides 

50  - 1000 

Convection  forcée 

Dans  les  gaz 

25  - 250 

Dans  les  liquides 

50  - 20000 

Con  vection  avec  chansement  de  phase 

Ebullition  ou  condensation 

2500  - 100000 

> Exercice  d’application  2: 

1.  De  l’air  froid  à Tm  = 16 °C,  soufflé  par  un  ventilateur,  s’écoule  sur  un  mur  de  10  m2  de 
surface,  maintenue  à une  température  constante  de  23 °C.  Si  le  coefficient  de  transfert  de 
chaleur  par  convection  est  de  100  W/m2.°C,  déterminer  la  puissance  thermique  échangée 
entre  l’air  et  le  mur. 

2.  L’air  à Tœ  = 30°C  est  en  contact  avec  une  plaque  plane,  de  2 m2  de  surface,  maintenue  à une 
température  constante  de  15°C.  Si  le  coefficient  de  transfert  de  chaleur  par  convection  est  de 
10  W/m2.°C,  déterminer  la  puissance  thermique  échangée  entre  le  fluide  et  la  plaque. 

3.  Du  pétrole,  à la  température  Tœ  = 120°C,  est  mis  en  circulation  forcée  dans  un  tube 
cylindrique  de  diamètre  D = 1 cm  et  de  température  Ts  = 60 °C.  Le  coefficient  de  transfert 
de  chaleur  par  convection  forcée  est  h = 1740  W/m2.°C.  Calculer  le  flux  de  chaleur 
échangé  entre  le  fluide  et  le  tube,  par  unité  de  longueur  du  tube. 

Solation: 

1.  Le  flux  de  chaleur  échangé  par  convection  forcée  entre  l’air  froid,  soufflé  par  le  ventilateur,  et 
le  mur  donné  par  la  loi  de  Newton  : 

O = h.S.  (T s - Tœ)  * 7 KW 


2.  Le  flux  de  chaleur  transféré  par  convection  de  la  plaque  à l’air  se  détermine  grâce  à la  loi  de 
Newton: 

O = h.  S.  (Tœ  - T g)  ~ 300  W 
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3.  Le  flux  de  chaleur  transféré  par  convection  du  pétrole  au  tube  se  calcule  en  appliquant  la  loi 
de  Newton  : 

O = h.S.  (Tœ  - T s) 


où  S est  l’aire  de  la  surface  latérale  intérieure  du  tube  (de  longueur  L ):  S = n.D.  L. 
Il  s’ensuit  que  le  flux  de  chaleur  échangé  par  unité  de  longueur  du  tube  est: 

O'  = | = h.  (Tœ  — T s)  * 3280  W/m 


> Exercice  d’application  3: 

Un  réfrigérateur  est  placé  dans  une  salle  où  la  température  de  l’air  est  de  21  °C.  Les  parois  de  ce 
réfrigérateur  sont  supposées  avoir  une  épaisseur  e = 3 0 mm  et  une  conductivité  thermique 
équivalente  de  0, 1 W/m.  °C.  La  température  sur  les  surfaces  externes  de  cet  appareil  est  égale  à 16°C 
et  le  coefficient  de  transfert  de  chaleur  à l’extérieur  est  de  10  W/m2.  °C.  Calculer  la  densité  du  flux  de 
chaleur  transféré  à travers  les  parois  de  ce  réfrigérateur.  En  déduire  la  température  sur  ses  surfaces 
internes. 

Solution: 

■ Paroi  du  réfrigérateur:  (À,  = 0, 1 W/m.  °C,  e = 30  cm) 
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La  densité  du  flux  de  chaleur  transféré  par  convection  de  l’air  vers  la  surface  frontière  externe  de  la 
paroi  du  réfrigérateur  est  donné  par: 


*Pcv  h X (T œ Tpe') 


En  régime  permanent,  la  densité  du  flux  de  chaleur  transféré  par  conduction  à travers  une  paroi 
simple  (milieu  conductif  sans  source  ni  puits  de  chaleur)  d’épaisseur  e et  dont  les  surfaces  frontières 
sont  portées  aux  températures  T1  et  T2,  s’exprime  comme  suit  : 


Bilan  thermique  sur  la  paroi: 

— *Pcv  — *Pcd 


*Pcd  A,  X 


T - T 

1 pe  ‘pi 


<p  = hx  (7*00  - Tpe)  = Ax 


T -T  ■ 

1 pe  1 pi 


Par  conséquent: 

■ La  densité  du  flux  de  chaleur  transféré  à travers  les  parois  du  réfrigérateur  est  : 

(p  = hx  (Tœ  — Tpe)  = 50  W /m2 

■ La  température  sur  les  surfaces  frontières  internes  Tpi  est  : 

<p=AxTpe~Tpi  =>  Tpi  = Tpe-jX<p  = 1°C 

3.3.  Rayonnement: 

Le  phénomène  de  rayonnement  thermique  constitue  un  mode  particulier  de  transfert  thermique  dans 
lequel  les  porteurs  d’énergie  ne  sont  plus  représentés  par  des  particules  de  substance,  mais  par  des 
ondes  électromagnétiques.  Le  mécanisme  de  rayonnement  est  créé  par  l’émission  et  l’absoiption  des 
ondes  électromagnétiques  porteuses  d’énergie  rayonnante  et  par  la  transformation  de  celle-ci  en 
chaleur  avec  le  changement  de  l’état  énergétique  (énergie  interne)  des  coips  qui  la  retiennent. 

Tous  les  corps  émettent  continûment  de  l’énergie  par  leur  surface,  sous  forme  d’un  rayonnement 
d’ondes  électromagnétiques  et  ce  d’autant  plus  que  leur  température  est  élevée.  Inversement,  soumis  à 
un  rayonnement  incident  ils  en  absorbent  une  partie  qui  se  transforme  en  chaleur  (stockée  en  leur 
sein).  L’intensité  du  rayonnement  émis  ou  absorbé  dépend  de  la  température  et  la  nature  du  matériau 
(propriétés  d’émission  et  d’absoiption  des  matériaux). 

Le  rayonnement  est  un  processus  physique  de  transmission  de  chaleur  entre  deux  corps  de 
températures  différentes,  lorsqu’ils  sont  séparés  par  un  milieu  matériel  transparent  ou  senti 
transparent  ou  même  lorsqu’un  vide  existe  entre  eux.  Donc,  contrairement  à la  conduction  et  à la 
convection,  ce  mode  de  transfert  de  chaleur  ne  nécessite  pas  la  présence  d’un  support  matériel  solide 
ou  fluide. 

> Exemples: 

► Rayonnement  reçu  par  les  parois  internes  d’une  chambre  de  combustion  ou  d’un  four 
industriel  en  provenance  des  gaz  chauds. 

► Rayonnement  échangé  entre  les  parois  d’un  four,  d’une  chaudière... 

> Remarque: 

Aux  températures  ordinaires , le  transfert  de  chaleur  par  rayonnement  est  négligeable  (la  conduction  et 
la  convection  jouent  un  rôle  important),  mais  il  peut  devenir  notable  et  prépondérant  lorsque  le  niveau 
de  température  augmente  (T  > 400 °C). 
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> Exercice  d’application  4: 

Une  plaque  horizontale  (A  = 30  W/m2.°C)  de  dimensions  900  mm  x 600  mm  x 30  mm  est 
maintenue  à Tps  = 300°C.  L’air  à Tœ  = 30°C  s’écoule  sur  la  paroi  supérieure  de  cette  plaque.  Si 
250  W de  chaleur  est  perdue  par  rayonnement  par  la  paroi  supérieure  de  la  plaque  et  le  coefficient  de 
transfert  de  chaleur  par  convection  air-paroi  supérieure  vaut  22  W/m 2 . °C,  quelle  sera  la  température, 
notée  Tp  i,  sur  la  paroi  inférieure  de  la  plaque. 

Solution: 

• Plaque:  (A  = 30  W/m.  °C , e = 30  mm,  l = 600  mm,  L = 900  mm,  TPi  =?  et  TPs  = 
300°C) 

■ Air  environnant  en  écoulement  forcé:  (Tœ  = 3 0 C et  h = 22  W/m2.  °C) 

■ Puissance  thermique  perdue,  par  rayonnement,  par  la  plaque  : 0rnet  = 250  W 


Bilan  thermique  sur  la  paroi: 

^cd  = &CV  "h  @r,net 


(9) 


Avec  : 

0cd  : le  flux  de  chaleur  transmis  par  conduction  à travers  la  plaque  ; 


&cv\  le  flux  de  chaleur  échangé  par  convection  entre  la  surface  frontière  supérieure  de  la  plaque  vers 
l’air  environnant; 


0 met '■  Ie  fh|x  de  chaleur  perdu  par  rayonnement  par  la  surface  frontière  supérieure  de  la  plaque. 
D’après  la  loi  de  Fourier,  on  a : 


0cd  = Ax  Sx 


T — T T — T 

1 P,l  1 V,S  „ r , 1 P, l 1 p, s 

— — = Ax  Lxlx  — — 


D’après  la  loi  de  Newton , on  a : 

0CV  = hx  S x ( TPiS  -Tœ)  = hx  Lxlx  ( TP:S  - Tœ ) 


L’équation  du  bilan  (9)  devient: 

T —T 

AxLxlx  ^ = h x L x l x (Tps  - Tœ ) + 


0, 


r,net 


D’où  on  tire  l’expression  de  la  température  Tp  i : 


WWVfc 
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Chapitre  2:  Rayonnement  Thermique 


1.  Introduction: 

Le  phénomène  de  rayonnement  thermique  constitue  un  mode  particulier  de  transfert  thermique  dans 
lequel  les  porteurs  d’énergie  ne  sont  plus  représentés  par  des  particules  de  substance,  mais  par  des 
ondes  électromagnétiques.  Ainsi,  contrairement  à la  conduction  et  à la  convection,  le  rayonnement  ne 
nécessite  pas  la  présence  d’un  support  matériel  solide  ou  fluide. 

2.  Généralités  et  définitions 

2.1.  Rayonnement  thermique: 

Notre  intérêt  est  porté  uniquement  sur  le  rayonnement  thermique  qui  est  celui  émis  par  la  matière  du 
fait  de  sa  température  et  au  détriment  de  son  énergie  calorifique.  11  s’étage  de  l’ultraviolet  à 
l’infrarouge:  [0, 1 — 100  ] fjm.  Cet  intervalle  spectral  se  découpe  ainsi:  0, 1 — 0, 4 fim  début  de 
l’ultraviolet,  0, 4 — 0,  8 pctn  domaine  visible,  0, 8 — 100  pcm  infrarouge. 

2.2.  Classification  des  corps  soumis  à un  rayonnement: 

Lorsqu’un  rayonnement  incident  (<P,  ) arrive  sur  la  surface  d’un  corps  semi  transparent,  il  est  en  partie 
réfléchi  (&r),  en  partie  absorbé  ( 0a ) et  le  reste  est  transmis  (Qt)  comme  l’indique  la  figure  ci-après. 


La  relation  entre  les  énergies  incidente,  réfléchie,  absorbée  et  transmise  s’écrit: 

0i  = 0r  + @a+  0t  ( Principe  de  conservation  de  V énergie) 

Elle  s’écrit  encore: 


p + a + t = 1 (1) 


ou: 


0r 


P — — est  le  coefficient  de  réflexion  ou  réflectivité , c’est-à-dire  la  fraction  du  rayonnement 


incident  réfléchi  par  le  corps; 


Oa 


Ct  — — est  le  coefficient  d’absorption  ou  absorptivité,  c’est-à-dire  la  fraction  du  rayonnement 


incident  absorbé  par  le  corps; 
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0t 

T — — est  le  coefficient  de  transmission  ou  transmitivité , c’est-à-dire  la  fraction  du  rayonnement 
incident  transmis  par  le  corps. 

Ces  paramètres  caractérisent  le  comportement  du  corps  vis-à-vis  du  rayonnement  qu’il  reçoit. 

On  distingue  trois  types  de  milieux  matériels  vis-à-vis  du  rayonnement:  deux  cas  limites  importants, 
les  milieux  parfaitement  transparents  et  les  milieux  opaques,  et  le  cas  général  des  milieux  semi 
transparents. 

Les  milieux  parfaitement  transparents  transmettent  tout  le  rayonnement  incident,  et  ce  quelle 
que  la  longueur  d’onde  X et  la  direction  de  ce  rayonnement.  Pour  ces  corps  la  relation  (1)  se 
réduit  à:  t = 1.  Le  vide  et  les  gaz  simples  (02,  N2,  H2  ...)  correspondent  à ce  type  de  milieu. 

Les  milieux  opaques  ne  transmettent  aucune  fraction  du  rayonnement  incident;  ce  dernier  est 
alors  partiellement  réfléchi  et  absorbé  (le  rayonnement  est  absorbé  sous  forme  de  chaleur  près 
de  l’impact).  Pour  ces  corps  la  relation  (1)  se  réduit  à:  p + a = 1.  La  majorité  des  liquides  et 
des  solides  sont  des  corps  opaques.  Ces  corps  arrêtent  la  propagation  de  tout  rayonnement  dès 
leur  surface  et  se  réchauffent,  de  ce  fait,  par  absorption  du  rayonnement. 

Les  milieux  semi  transparents  ou  partiellement  transparents  sont  caractérisés  par:  p + a + 
t = 1.  Dans  ce  cas,  la  propagation  du  rayonnement  s’accompagne  d’une  diminution  de 
l’énergie  transportée,  ce  qui  accroît  d’autant  l’énergie  interne  du  milieu  traversé.  C’est  le  cas 
de  certains  gaz  composés  ( C02 , H20,  CO,  CH4  ...)  et  de  certains  solides  (plastiques, 
verres...). 

Les  grandeurs  radiatives  p,  a et  t dépendent  non  seulement  de  la  nature  du  matériau,  de  son  épaisseur 
et  de  son  état  de  surface,  mais  également  de  la  longueur  d’onde  X et  de  la  direction  du  rayonnement. 

2.3.  Classification  des  grandeurs 
Les  grandeurs  radiatives  sont  distinguées  selon  : 

r La  composition  spectrale  du  rayonnement 

Si  la  grandeur  correspond  à l’ensemble  du  spectre  du  rayonnement  thermique,  elle  est  dite 

totale. 

Si  elle  concerne  un  intervalle  spectral  étroit  dX  centré  autour  d’une  longueur  d’onde  X,  elle 
est  dite  monochromatique . On  l’affecte  d’un  indice  /(pour  la  distinguer  de  la  grandeur  totale: 

gâ 

y La  distribution  spatiale  du  rayonnement 

Si  la  grandeur  correspond  à l’ensemble  des  directions  de  l’espace  dans  lesquelles  un  élément 
de  surface  peut  rayonner  ou  recevoir  du  rayonnement,  elle  est  dite  hémisphérique. 

Si  elle  caractérise  une  direction  donnée  de  propagation  du  rayonnement,  relativement  à la 
surface  considérée,  elle  est  dite  directionnelle . On  l’affecte  d’un  indice  ox  pour  la  distinguer 
de  la  grandeur  hémisphérique:  Gox 

3.  Corps  noir  et  lois  de  rayonnement: 

3.1.  Corps  noir 

Un  corps  qui  absorbe  tout  le  rayonnement  qu’il  reçoit,  quelle  que  soit  sa  longueur  d’onde  et  sa 
direction,  est  appelé  corps  noir.  Il  sert  de  référence  à partir  de  laquelle  les  propriétés  radiatives 
(d’émission  et  de  réception  du  rayonnement)  des  autres  corps  se  situent  en  retrait.  Son  coefficient 
d’absorption  «est  égal  à 1 et  ses  coefficients  de  réflexion  pet  de  transmission  rsont,  par  conséquent, 
nuis. 
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Le  corps  noir  est  aussi  l’émetteur  idéal,  il  a la  propriété  d’émettre,  à chaque  température  et  longueur 
d’onde,  le  maximum  d’énergie  thermique  pouvant  être  rayonné  par  la  matière. 


3.2.  Lois  du  rayonnement: 

Les  lois  du  rayonnement  que  nous  allons  établir  décrivent  l’émission  du  coips  noir.  Dans  tout  ce  qui 
suit,  les  grandeurs  relatives  au  corps  noir  seront  affectées  de  l’exposant”0”  pour  les  distinguer  de 
celles  relatives  au  corps  réels. 


3.2.1.  Loi  de  Planck:  émittance  monochromatique  du  corps  noir 
La  puissance  thermique  monochromatique  rayonnée,  par  unité  de  surface,  par  un  corps  noir  à la 
température  T et  à la  longueur  d’onde  A est  donnée  par  la  loi  de  Planck : 


M°À(A,T)  = 


(W/  2 ) (2) 


A5  x (exp Q x2 l) 

Les  valeurs  des  constantes  physiques  Ct  et  C2  sont  données,  en  unités  S.I.,  dans  le  tableau  suivant: 


T 

A 

Cx  = 2nx  h x C2 

hx  C 
C>=  K 

K 

(m) 

3,741  x 1(T16  W.m2 

0,014388  m.  °K 

W/m3 

Kelvin  ( °K ) 

(ptn) 

3,741  x 108  W.pm4/m2 

14388  pm.°K 

W/m2.pm 

La  fonction  A/8,  appelée  émittance  monochromatique  du  corps  noir , est  représentée  par  le  réseau  de 
courbes  ci-après. 
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Figure  2:  Distribution  spectrale  de  l’émittance  du  corps  noir  en 
fonction  de  la  température  absolue  (Loi  de  Planck). 


On  constate  qu’à  chaque  température  T correspond  une  courbe  ayant  un  maximum  situé  à une  valeur 
À,nax  de  la  longueur  d’onde  variable  avec  T.  Ces  courbes  présentent  une  dissymétrie  très  prononcée: 
leur  croissance  avec  A (courtes  longueur  d’onde)  est  beaucoup  plus  rapide  que  leur  décroissance.  Pour 
A < 0,  5 ■Amax  il  n’y  a pratiquement  plus  d’énergie  rayonnée  (moins  de  1%)  alors  qu’il  faut  atteindre 
A < 5.  Amax  pour  obtenir  le  même  résultat  dans  l’infrarouge. 

Ainsi,  dans  la  pratique  on  peut  considérer  que  l’étendue  utile  du  spectre  est  comprise,  à chaque 
température  T , entre  environ  0,  5 X A,nax  et  5 X A,nax. 

3.2.2.  Lois  de  Wien 

Ces  deux  lois  fournissent  respectivement  l’abscisse  et  l’ordonnée  du  maximum  d’émittance 
monochromatique  du  corps  noir  pour  chaque  température. 

3.2.2. 1.  Première  loi  de  Wien  ou  loi  de  déplacement 

Cette  loi  permet  d’exprimer  commodément  que  les  maxima  Amux  des  coui'bes  de  (fig.  2)  se 
déplacent  vers  les  petites  longueurs  d’onde  quand  la  température  augmente,  selon  la  relation: 

^max  X T = 2898  (/ um.  K ) (3) 

3.2. 2.2.  Deuxième  loi  de  Wien 

Cette  loi  fournit  l’émittance  monochromatique  maximale  M^max: 
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Mtawx  = b XTS  (xym2  (lm)  (4) 
avec  B une  constante  qui  a pour  valeur:  B = 1,  287  X 10-11  (W/m2.  pm.  K 5) 


3.2.3.  Loi  de  Stefan-Boltzmann  : 

La  puissance  totale  rayonnée,  par  unité  de  surface,  par  un  corps  noir  à la  température  T est  fournie  par 
la  loi  de  Stefan- Boltzmann  : 

M°(T)  = o x T4  ( W/m 2)  (5) 

, „ , „ , f <7  = 5,67xl(T8  W/m2. K4 

ou  a est  la  constante  de  Stefan- Boltzmann:  { „ .. 

(a  = 4, 88  x 1(T8  Kcal/h.m2.K 4 

M°(T ) est  appelée  émittance  totale  du  rayonnement  du  corps  noir. 


3.2.4.  Fraction  de  l’émittance  totale  contenue  dans  un  intervalle  spectral  donné:  Fa1t-a2t 

La  fraction  de  l’énergie  émise  par  une  surface  d’un  coips  noir  de  température  T dans  l’intervalle 
spectral  [A\,  A2]  est  donnée  par: 


F AxT-faT 


M°À(A,T)  x dA 
J0+ooM»a,r)  x dA 


(6) 


que  l’on  peut  également  écrire: 

1 

F axt-a2t  — g 


(A,  T)  xdA- 


rAl 

I M°À(A,  T ) x dA 
o 


(7) 


Soit  encore: 


Fa-iT-AïT  ~ Fq-a2t  ~ Fq-^t  (8) 


Pour  pouvoir  calculer  cette  fraction,  il  serait  intéressant  de  disposer  d’une  table  où  F 0_xr  est  donnée 
en  fonction  de  A.  T (voir  annexe  1). 


> Exemple  d’application  1: 

Le  soleil  rayonne  comme  un  corps  noir  à T=5800  K. 

1.  Calculer  l’émittance  totale  du  soleil,  M°. 

2.  Calculer  l’émittance  monochromatique  maximale  M*lmax  et  la  longueur  d’onde 
correspondante  A,mix. 

3.  Evaluer  la  fraction  de  l’émittance  totale  du  soleil  F située  dans  le  domaine  visible 

[0,4-0, 8]  pan. 

4.  On  a situé  l’étendue  utile  du  spectre  rayonné  par  le  corps  noir  à une  température  T entre  les 
longueurs  d’onde  0,5xAnmx  et  5xA,m,x.  A quel  pourcentage  de  l’énergie  totale  rayonnée  par  le 
corps  noir  à la  température  T cela  correspond-il? 

Solution: 

1.  Emittance  totale  du  soleil:  M°(T)  = a.  T4  « 64164  KW/m 2 

2-  et  A,,mx  sont  données  par  les  deux  lois  de  Wien: 

( 2898 

k max  = -Zjr-  ~ 0,  5 pm  infra  - rouge 
( = » 84473  (KM7m2.ftm) 
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3. 


On  calcule  d’abord  les  paramètres  d’entrée: 

ai  t = 
U2t  = 


2320  j um.  K 
4640  j um.  K 


La  table  fournit  les  valeurs  suivantes: 

f Fo-âxT  - 
(Fo-a2t 


0,  12  1 

0,5851 


=>  ^AiT-AzT  — 0,4611 


Environ  46%  de  l’énergie  solaire  sont  rayonnés  dans  le  visible  [0,4-0, 8 ] /Mil. 

4.  Calculons  F0,5xUT-5xwr: 

Les  paramètres  d’entrée  de  la  table  fournissant  les  Fü_xr  sont: 

«■  5 X K,aJ  = 1450  !«"■  K jWr  = 0, 0096 
5 x A,,,,,,/  = 14500  nm.K  ( Fo-s** ,,  = 0, 9666  o.sx 

C/C:  L’étendue  utile  du  spectre  rayonné  par  le  soleil  correspond  donc  à 95,7%  de  l’énergie  totale 
rayonnée  par  le  corps  noir  à la  température  considérée. 


4.  Rayonnement  des  corps  réels: 

4.1.  Emission  des  corps  réels 

Les  propriétés  radiatives  d’émission  des  coips  réels  sont  définies  par  rapport  à celles  du  corps  noir 
dans  les  mêmes  conditions  de  température  et  de  longueur  d’onde.  Ainsi,  on  définit  l’émissivité  f d’un 
corps  réel  comme  le  rapport  de  son  émittance  à celle  du  corps  noir. 

Les  coefficients  d’émission  ou  émissivités  monochromatique  ex  et  total  s d’un  coips  réel  sont  définis, 
respectivement,  par: 


ex  a,T)  = 


M°U  70 


(9  — a) 


et 


M(T ) 


(9 -b) 


MX(A,  T ) et  M(T)  sont,  respectivement,  les  émittances  monochromatique  et  totale  du  corps  réel. 


L’émissivité  d’un  corps  réel  peut  dépendre  également  de  la  direction  ; elle  est  alors  appelée  émissivité 
directionnelle  et  notée  eox. 


> Exemple  d’application  2:  (Extrait  de  TD  de  GP-IGM2_  2012/2013) 

Déterminer  l’émissivité  totale  hémisphérique  £ et  l’émittance  totale  M d’une  surface  réelle  rayonnant 
à T = 838  °C  dont  l’émissivité  hémisphérique  monochromatique  £x  est  de  0,8  pour  les  longueurs 
d’onde  inférieures  à 1,  5 pm,  de  0,  6 pour  les  longueurs  d’onde  de  1,  5 à 2,  5 j um  et  de  0, 4 pour  les 
longueurs  d’onde  supérieures. 


Solution: 

On  considère  une  surface  rayonnant  à T = 838°C  et  dont  l’émissivité  hémisphérique 
monochromatique  £x  est  constante  sur  trois  bandes  spectrales  : 

(£x  = 0, 8 (0  < A < 1,5  (j um)) 

£2  = 0, 6 (l,  5 < A < 2,5  (j um)) 

£3  = 0,4  ( A>  2,5  (j um)) 

Une  telle  surface  est  appelée  surface  semi-grise  ou  surface  grise  par  bandes. 
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Calculons  l’émissivité  totale  de  cette  surface  pour  T = 838°C:  e(T  = 838°C) 


On  rappelle  que  l’émissivité  monochromatique  £^(Â,  T ) et  l’émissivité  totale  e(T)  sont  définies 
respectivement  par  : 

E^À’T^>  -Ta, 


et 


£(J)  = 


M'iUT) 
M(T ) 


M°(T) 

► L’émissivité  totale,  à la  température  T , de  la  surface  rayonnante  (corps  réel)  est  donnée  par  : 

M{T) 


M(T)  1 f+0°  1 f+0° 

£(r)  = MW  = Mwxj0  "MndÀ=üsm*l 


On  pose  : 


Donc  : 


A1  = 1,  5 j um  et  Â2  = 2,5  fim 


e(T)  = 


axT4 


x 


rA±  rA 2 r+OO 

£±x  M%UT).dA  + £ 2x\  M°À(A,T).dA+£3x  M°À(A,T).dA 

Jo  Ja,  h2 


D’où  : 


£(T)  — £i  x F o-Al  T + £2  X F A1.T-A2.r  + £3  X (l  — Fo-A2.t) 
Déterminons  les  valeurs  de  F o-^.r  et  F 0_^2  j 

Les  paramètres  d’entrée  de  la  table  donnant  ces  fractions  sont  : 

f/li  x T = 1666,  5 j um.  K ~ 1670  j um.  K 
U2xT  = 2777,  5 ptm.  K » 2780  / um.  K 


D’après  la  table,  on  a: 

F o_Z7(A  T = 1660)  + F0_ÀT(Z  T = 1680)  0,  0247  + 0,  0266 


11  s’en  suit  que  : 


°-^T  ~ 2 ~~  2 

Fo -a2.t  = 0,  2234 

£{T  = 838°C)  = 0,45 

► Emittance  totale  de  surface  rayonnante  : 

M(T ) = £ (T)  x M°(r)  = £{T ) x a x T4 

A.N.  : M(T  = 838°C)  = 38873  W * 39  KW 


= 0,02565 


4.2.  Corps  gris: 

Un  corps  gris  est  un  corps  réel  dont  les  propriétés  radiatives  (d’émission  et  de  réception)  sont 
indépendantes  de  la  longueur  d’onde  A,  à une  température  donnée: 

(£A,ox  £ox  (aA,ox  ~ > aox 

1 £a  -»  £ l aÀ->  a ■ 

Le  corps  gris  est  un  cas  particulier  (usuel)  comme  les  cas  suivants: 
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■ Corps  réel  et  diffusant  en  émission  et  en  réception  ( isotrope ) : émissivité  et  absorptivité 
indépendantes  de  la  direction  de  propagation  du  rayonnement. 

f£Â,ox  £Â  (aA,ox  ax 

1 £ox  -»  £ l Uox  -»  « 

■ Corps  gris  et  diffusant  en  émission  et  en  réception:  émissivité  indépendante  de  la  longueur 
d’onde  et  de  la  direction  de  propagation  du  rayonnement. 

£Â,ox  £ et  aÂ,ox  a 


4.3.  Loi  de  Kirchhoff: 

La  loi  de  Kirchhoff  indique  que  les  coefficients  d’émission  et  d’absoiption  sont  égaux  pour  la  même 
longueur  d’onde  et  la  même  direction  de  propagation  du  rayonnement  émis  par  la  surface  d’un  corps 
ou  incident  sur  celle-ci. 

£À,ox  tX ^ox  (10) 

Les  corps  réels  peuvent,  en  général,  être  assimilés  à des  corps  gris  et  diffusant  (propriétés  radiatives 
indépendantes  de  la  direction).  Un  corps  réel  est  alors  caractérisé  par  une  seule  grandeur  radiative: 
son  émissivité  £. 

a = £ (11) 

Relations  générales  entre  U absorptivité  et  l’ émissivité 
■ Pour  un  corps  réel  quelconque:  £ 4=  a 

• Pour  un  corps  gris  et  diffusant:  £%  = a^  = £ = a < 1 

• Pour  un  corps  noir:  £ = a = 1 

Dans  le  tableau  suivant,  on  donne  l’émissivité  totale  hémisphérique  £ pour  quelques  corps  (à  T = 
20  °C). 


Matériaux 

Etat  de  surface 

£ 

Aluminium 

Poli 

0,06 

oxydé 

0,30 

Cuivre 

Poli 

0,04 

Très  oxydé 

0,75 

Or 

Pur  très  poli 

0,02 

Argent 

poli 

0,02 

Zinc 

Pur  très  poli 

0,02 

Galvanisé 

0,2-0, 3 

Brique 

Ordinaire 

0,93 

Réfractaire  (à  1000° C) 

0,6 

Carbone  ( dépôt  de  noir  de 
fumée) 

0,95-0,98 

Verre  ordinaire 

0,94 

Marbre 

0,95 

Peintures 

Aluminium 

0,30-0,60 

Noire 

0,98 

Blanche 

0,90 

Peintures  à l’huile  toutes 
couleurs 

~0, 90 

Papier 

0,95 

www . iesuisetudiant . coin 
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5.  Echanges  radiatifs  entre  surfaces  séparées  par  un  milieu  parfaitement  transparent: 

On  rappelle  qu’un  milieu  matériel  parfaitement  transparent  (r  = 1)  est  un  milieu  non  absorbant  (ni 
émissif.  a = e = 0)  et  ne  participe  donc  pas  aux  échanges  radiatifs.  L’air  et  les  gaz  simples  en  sont 
des  exemples. 

On  se  propose  d’étudier  les  échanges  radiatifs  entre  surfaces  séparées  par  un  milieu  parfaitement 
transparent.  Cette  situation  est  très  largement  rencontrée  dans  la  pratique.  En  effet,  on  trouve  de 
nombreuses  applications  industrielles  pour  lesquelles  les  échanges  radiatifs  entre  surfaces  sont  peu 
affectés  par  le  milieu  transparent  qui  les  sépare. 

Les  surfaces  seront  d’abord  supposées  noires,  ensuite  grises  et  diffusantes  en  émission  et  en 
réflexion. 

Hypothèses  simplificatrices  : 

Toutes  les  surfaces  considérées  seront  supposées  opaques,  isothermes,  homogènes  et  diffusantes. 


5.1.  Echanges  radiatifs  entre  surfaces  noires: 

5.1.1.  Facteur  déformé  géométrique  entre  deux  surfaces  F y: 

Considérons  deux  surfaces  noires  S j et  Sj  maintenues  à des  températures  uniformes  T,  et  T y et 
pouvant  rayonner  mutuellement  l’une  vers  l’autre.  Soit  0t  = Sj  X M-*  le  flux  total  émis  par  Sj,  dont 
seule  une  partie  <Z>,y  est  reçue  par  Sy. 


On  appelle  facteur  déformé  de  Sj  vers  Sy  le  rapport: 


Ainsi: 


Fa  = 


0i 


(12) 


0 < Fij  < 1 


Figure  3:  Echanges  radiatifs  entre  deux  surfaces  Stet  Sj. 


F ij  est  un  nombre  sans  dimension  représentant  la  fraction  du  rayonnement  émis  par  la  surface  Sj  dans 
toutes  les  directions  interceptées  par  la  surface  Sj.  Ce  facteur  de  forme  ne  dépend  que  de  la  géométrie 
et  de  la  disposition  relative  des  surfaces  rayonnantes. 


5.1.2.  Propriétés  des  facteurs  déformé 
5. 1.2.1.  Règle  de  réciprocité 

Les  facteurs  de  forme  entre  deux  surfaces  Sj  et  Sy  vérifient  la  relation  suivante  désignée  sous  le  nom 
de  règle  de  réciprocité : 


WWVf 


Sj'  F îj  = Sy.  Fji 


(13) 


nement  Thermiqu 
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5. 1.2.2.  Règle  de  complémentarité  ou  d’addition: 

Considérons  une  enceinte  fermée  constituée  de  N surfaces  noires  Sx,  S2,  ....  Si,...  SN, 
individuellement  isothermes  (figure  4).  La  définition  du  facteur  de  forme  géométrique  induit 
immédiatement  la  règle  de  complémentarité'. 

On  a: 

N N N 

^ ®ij  = x^)  = ^x^  Fij  avec  1 = 1,2.. .N 

;=i  7=i  7=i 


Donc: 


N 

N 

N 

N 

II 

XI 

II 

c^r 

X 

• = Î/'7  = - 

- = YJFNj 

7=1 

7=1 

7=1 

7=1 

= 1 


C’est-à-dire,  la  somme  des  facteurs  de  forme  associés  à une  surface  S j (avec  i = 1,2  ...AO  est  telle 
que  : 

N 

^ F a = i (14) 

7=1 


Figure  4:  Enceinte  fermée  constituée  de  N surfaces  noires. 

5. 1.2.3.  Règle  de  décomposition 

Considérons  deux  surfaces  S)  et  S).  Si  S,  se  décompose  en  deux  surfaces  Si  a et  Sib  (figure  5),  on 
peut  écrire: 

Si- Fij  = Sj  a.  Fi  a j + Si  b.  Fi  b j (15) 


Figure  5: 

Ces  règles  sont  utiles  pour  la  détermination  des  facteurs  déformé  géométrique  de  plusieurs  surfaces 
en  présence. 


> Remarque: 

■ Le  facteur  de  forme  Fa  d’une  surface  S,  sur  elle-même  est  nul  lorsque  cette  surface  est 
plane  ou  convexe. 
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■ F^e st  différent  de  0 quand  S j est  concave. 

5.1.3.  Quelques  facteurs  déformés  géométriques: 

Il  existe  dans  la  littérature  des  tables  et  des  abaques  qui  fournissent  des  valeurs  approchées  des 
facteurs  de  forme  géométriques  Fy.  On  donne  ici  quelques  formules  utiles  permettant  de  calculer  F y 
dans  de  le  cas  de  configurations  simples. 

► une  surface  convexe  S2  à l’intérieur  d’une  cavité  (figure  6): 


► Deux  surfaces  délimitant  un  volume  clos,  l’une  n’étant  pas  concave  5/  (figure  7): 


► Plans  parallèles  proches  (figure  8):  effets  de  bord  négligeables 


► Cylindres  parallèles,  emboîtés,  de  grande  longueur,  quelle  que  soit  la  distance  S de  leurs  axes 
(figure  9)  : 

R-t  R-t 

F n = 0,  F12  = 1,  F2i  = — — et  F22  = 1 — -i 

► Sphères  emboîtées,  quelle  que  soit  la  distance  S de  leurs  centres  (figure  9): 
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> Exemple  d’application:  (Voir  TD,  série  n°l) 

5.1.4.  Flux  radiatif  net  échangé  entre  surfaces  noires 

a)  Cas  de  deux  surfaces  noires  Sj  et  Sj 

Les  deux  surfaces  en  vis-à-vis  étant  noires,  donc  tout  rayonnement  reçu  par  l’une  ou  l’autre  est 
entièrement  absorbé  (il  n’y  a pas  de  réflexion  ni  de  transmission). 

Le  flux  radiatif  net  transféré  de  Sj  vers  Sj,  noté  Ornet.  j_>j,  s’obtient  en  faisant  un  bilan  entre  le  flux 
émis  par  Sj  vers  Sj  (perte  d’énergie  thermique  pour  Sj)  et  le  flux  absorbé  par  Sj  en  provenance  de  Sj 
(gain  d’énergie  thermique  pour  Sj).  Il  est  donné  par  la  relation: 

&r,net.  t-J  = % ~ Ojj  = Fij  X ~ Fji  X = $1-  Fij  X K ~ Sj.  F jt  X M°j 

En  tenant  en  compte  de  la  règle  de  réciprocité  et  de  la  loi  de  Stefan-Boltzmann,  l’expression  de 
Omet.  i->j  devient  : 

Omet.  i->j  = Sj.  Fjj  x « - M°j)  = Sj.  Fjj  x a x (7?  - Tf ) (16) 

On  remarque  que  : 

Or,net:  j->i  ~ ~ Or  nei . i_>y 

b)  Cas  d ’une  enceinte  noire  fermée 

Considérons  une  enceinte  fermée  constituée  de  N surfaces  noires  Si,  S 2.  ...,  Sj, ...  SN, 

individuellement  isothermes,  aux  températures  T i,  T 2,  ....  Tj,...  Tn  ( figure  4).  Dans  la  pratique,  on 
est  amené  dans  certaines  situations  à combler  les  espaces  vides  entre  les  surfaces  (constituant 
l’enceinte)  à l’aide  de  surfaces  fictives  isothermes  de  température  égale  à celle  du  milieu  ambiant  où 
se  trouve  placée  l’enceinte.  C’est,  par  exemple,  le  cas  d’un  four  fonctionnant  à porte  ouverte  ( voir 
TD). 

Nôtre  objectif  est  de  trouver  l’expression  du  flux  radiatif  net  échangé  par  chaque  surface  de  l’enceinte. 

Pour  cela,  faisons  un  bilan  des  échanges  radiatifs  sur  la  surface  Sj,  cette  surface  émet  un  flux  total 
et  absorbe,  en  provenance  de  toutes  les  surfaces  composant  l’enceinte,  y compris  elle-même  si  elle  est 
concave  (F a * 0),  des  flux  Ojj  (j  = 1,2  ...  N).  Le  flux  radiatif  net  échangé  par  la  surface  est: 

N N N 

Or, net, i =<**“£  Ojj  = Ot  ~ ^ Fjj.  ®j  = Sj.  M°  - ^ Fjj.  Sj.  M°j 

7=1  7=1  7=1 

soit,  compte  tenu  des  règles  de  réciprocité  et  de  complémentarité: 

(N  \ N N 

M°j  ~ £ F jj.  M°j  = Sj  X £ Fjj.  (M°  -M°j)  = Yj  Omet.  t-*j  (17  - a) 

7=i  B J 7=1  « i 

www. lesuisetudiant . com 

Cours  de  Génie  des  Procédés-Transfert  de  Chaleur  Chapitre  2:  Rayonnement  Thermique 

12 


http://www.jesuisetudiant.com 


Un  assistant  Scolaire  Polyvalent 


http://www.jesuisetudiant.com  N'hesitez  Pas  de  Visiter  Notre  Site  Un  assistant  Scolaire  Polyvalent 

FST-Marrakech-Dévartement  de  Physique  Appliquée F.I.:  Génie  des  Matériaux-C.1.2  f 2014/2015 1 


soit,  compte  tenu  de  la  loi  de  Stefan- Boltzmann  : 

N N 


0 


•,net,i  = ^ 0r.net.  i->j  = St  X <T  X £ F y X (7?  - T) ) (17  - b ) 

7=1  7=1 


5.1.5.  Représentation  des  échanges  radiatifs  entre  surfaces  noires  par  analogie  électrique 
On  remarque  que  la  relation  exprimant  le  flux  échangé  entre  deux  surfaces  noires  Si  et  Sj 


0 


r,net\  i-*j 


i = St.  F ij  x (Mf  — Mj  ) = 


M°  - M°  Ml  - M°j 


(18) 


Si- F a 


Sj-  F ji 


est  analogue  à la  relation  exprimant  le  courant  électrique  /y  qui  circule  dans  un  réseau  électrique 
entre  deux  nœuds  de  potentiels  Vt  et  V j séparés  par  une  résistance  Ry. 

V - V 

v i v j 


hi  Ru 


L’échange  radiatif  entre  deux  surfaces  noires  5,  et  Sj  peut  donc  être  représenté  à l’aide  du  schéma 
électrique  analogique  suivant  (figure  10): 


M°  Sl.Flj  jyjO 

HI — Figure  10  : 


0 i 


'r.net:  l->2  1 

OU 


Sj-  F ji 


L’échange  radiatif  entre  trois  surfaces  noires  5j,  S2  et  53  constituant  une  enceinte  fermée  peut  être 
représenté  par  le  schéma  suivant  (figure  11): 

1 


■ Le  flux  total  émis  par  est  donné  par: 

3 

0r,net,l  I*  ,net\  l—>j  0r,net\  1->1  F 0r,net\  1-^2  "f  0r,net\  l->3  0r,net\  l->2  0r,net\  l->3 

7=1  =0 

0r,net,i  = $1-  F 12.  (M°  - Mf)  + Sv  F13.  (M?  - M°3)  = Si.  a.  (f12.  ( Ti  - Tj)  + F13.  (T*  - Tj)) 


■ Le  flux  total  émis  par  S2  est  donné  par: 
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3 

^ r,net,2  ^ ®r,net\  2 ~*j  ® r,net\  2->l  ® r,net\  2->2  "b  ®r,net.  2->3  ^ r,net : 2-»l  ~b  ®r,net\  2->3 

7=1  =0 

®r,net,2  = S2.F21.  (M°2  - M°)  + S2.F23.  (M°2  - M°3)  = S2.(J.  (f21.  (r|  - T\)  + F23 . (Tt  - T43)) 

■ Le  flux  total  émis  par  S3  est  donné  par: 

3 

®r,net,3  ^ ®r,net\  3 ->/'  ® r,net\  3->l  "b  ®r,net.  3->2  ~b  ®r,net\  3->3  ^ r,net : 3->l  ~b  ®r,net\  3-»2 

7=1  =0 

®r,net,3  = S3.F31.  (M3  - M°)  + S3.F32.  (M3  - M°2)  = S3.c T.  (f31.  ( T$  - T\)  + F32 . (Tl  - Tl)) 

> Exemples  d’application:  ( voir  TD,  série  n°l) 

5.2.  Echanges  radiatifs  entre  surfaces  grises: 

On  s’intéresse  à l’étude  des  échanges  radiatifs  entre  surfaces  grises  et  diffusantes  en  émission  et  en 
réflexion,  séparées  par  un  milieu  parfaitement  transparent  à toute  longueur  d’onde. 

5.2.1.  Notion  de  radiosité  d’une  surface  grise  (réelle) 

On  rappelle  que  si  est  le  flux  incident  sur  une  surface  opaque  d’aire  S,  seule  une  fraction  0a  est 
absorbée,  une  fraction  Qr  étant  réfléchie. 

h’ éclairement,  noté  E ,de  la  surface  S correspond  à la  quantité  : 

e = Y (W/m2)  (19) 

L ’émittance  de  cette  surface  est  : 

M = £XM°  = £Xffxr4  (w  / m 2) 

La  radiosité,  notée  /,  de  la  surface  5 est  le  flux  total  émis  par  cette  surface  (figure  12):.  Ce  flux  total 
résulte  pour  une  part  de  l’émission  propre  de  la  surface  (émittance  M)  et  pour  une  autre  part  de  la 
fraction  de  l’éclairement  qu’elle  réfléchit  (p  X E): 

J=M+pxE=£XM°+pxE  (W  /m2)  (20) 

Partie  réfléchie  de  ^ _ Radiosité  (J) 


Figure  12  : Partie  absorbée  de 

~ l’éclairement  (a.  E ) 


La  surface  en  question  est  opaque,  grise  et  diffusante  en  émission  et  réflexion,  donc: 

t = 0 et  p = 1 — a = 1 — £ 

Par  suite  : 

J = £ x M°  + (1  - £)  x E (W/m2)  b(21) 

www . lesuisetudiant . coin 

Cours  de  Génie  des  Procédés-Transfert  de  Chaleur  Chapitre  2:  Rayonnement  Thermique 

14 


http://www.jesuisetudiant.com 


Un  assistant  Scolaire  Polyvalent 


http://www.jesuisetudiant.com  N'hesitez  Pas  de  Visiter  Notre  Site  Un  assistant  Scolaire  Polyvalent 

FST-Marrakech-Dévartement  de  Physique  Appliquée F.I.:  Génie  des  Matériaux-C.1.2  f 2014/2015 1 


Le  flux  total  issu  d’une  surface  grise  diffusante  d’aire  S est  égal  à : 

0 = s x / = S x (f  x M°  + (1  - £)  x E)  ( W ) 

Le  flux  radiatif  net  Or  net  perdu  par  la  surface  grise  S est  égal  à la  différence  entre  le  flux  émis  et  le 
flux  absorbé  par  cette  surface,  soit  : 

Qrnet  = (S  x £ x M°)  - (5  x a X E)  = S X £ x (M°  - E)  ( W ) (22) 

Eliminons  l’éclairement  £ de  cette  relation  en  la  remplaçant  par  son  expression  tirée  de  (21),  on 
obtient  : 

( n J ~ £X  5 X £ . 

Or, net  = S X S X IM0  - ±_£  ) = ^ X (M0  - /)  (23) 


En  éliminant  M°  de  cette  relation,  on  obtient  : 

Or,net  = S X (J  - E)  (24) 

Ainsi,  le  flux  radiatif  net  perdu,  par  unité  de  surface,  par  une  surface  grise  S est  aussi  égal  à la 
différence  entre  sa  radiosité  J (densité  de  flux  de  chaleur  quittant  la  surface  par  émission  et  réflexion) 
et  son  éclairement  E (densité  de  flux  arrivant  sur  la  surface). 


On  dispose  donc  de  Pois  relations  distinctes  fournissant  le  flux  net  rayonné  par  une  surface  opaque 
grise  et  diffusante  S f 

Si  x Si 


O, 


,net.i  = Si  x Si  x (. M°i  - Ei ) = -i— - x{M°i-Ji)=SiX(Ji-  Ei ) (2  5) 

A ^ i 


5.2.2.  Représentation  des  échanges  radiatifs  entre  surfaces  grises  par  analogie  électrique 
L’expression  du  flux  net  rayonné  par  une  surface  grise  et  diffusante  5) 

Si  x Si  /.mr{l  w N M°  — Ji 


O , 


r.net.i 


1 -£, 


X 


(<-/t)  = ^T^i  (26) 


Si  x Si 


peut  être  interprétée,  en  termes  d’analogie  électrique,  à l’aide  d’un  circuit  dans  lequel  une  résistance 
ç — L , soumise  à une  différence  de  potentiel  — Ji,  est  traversée  par  un  courant  égal  à Orneti 
(figure  13). 


1 ~st 


M°i  si  X £i  f 

^ m lr— 


O, 


r.net.i 


Figure  13  : 


Considérons  deux  surfaces  grises  et  diffusantes  Sj  et  Sj  rayonnant  l’une  vers  l’autre.  Le  flux  total 
quittant  la  surface  5)  est  St  X Seule  la  fraction  de  ce  flux  atteint  la  surface  Sj,  donc  le  flux  issu 
de  Sj  et  arrivant  sur  Sj  est  Sj  X F jj  X Réciproquement,  le  flux  quittant  S) et  atteignant  Sj  est 
Sj  x F jj  X J j.  Le  flux  radiatif  net  échangé  entre  Sj  et  Sj  a pour  expression  : 

Or.net.i^j  =SjX  F jj  X {] \ - Jj)  (27) 


Cet  échange  peut  être  représenté  à l’aide  du  schéma  électrique  suivant  : 


Ji 


Sj  X F a 


9 ► 


WV5  ^ 


ni 
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Par  conséquent,  le  schéma  électrique  complet  représentant  les  échanges  radiatifs  entre  deux  surfaces 
grises  et  diffusantes  en  vis-à-vis  est  illustré  sur  la  figure  15. 


r,net,i 


- r,net : i—>j 


Si  X Ei 


Si-  Fij 


1 £j 


Sj  X Ej 


Figure  15: 


Dans  ce  cas,  on  a : 


M°  - M° 


SiXax(Tf-  Tj ) 


v r,net : i^>j 


i-£, . i , i-£,  i_, 

S-^  + SiTij+S^Ij  e,  1 + F,i+Slx{sl  1) 


5.2.2. 1.  Facteur  de  forme  gris  fq: 

Par  analogie  avec  l’échange  radiatif  entre  deux  surfaces  noires,  on  introduit  le  facteur  de  forme  gris, 
noté  'fij,  défini  par: 

r.  = —, \ (29) 


1 1 s 1 

--7  + — + ^ — -1 

s.  F..  S.  s. 
i y j \ j 


5.2.2.2.  Quelques  facteurs  de  formes  gris: 

Plans  parallèles  proches  (effets  de  bord  négligeables):  F 12  = 1 et  S1æS2- 

^«=1—7 — (30) 

— + 1 


S1  S2 


Cylindres  parallèles,  emboîtés,  de  grande  longueur,  quelle  que  soit  la  distance  S de  leurs  axes: 


5,  R, 


/ * » t 

(F]2  =7  et  — = — ) 

12  C L> 


*2 


*2  *2  v £ 2 


S.  R. 

Sphères  emboîtées,  quelle  que  soit  la  distance  <5 de  leurs  centres:  ( F n — 1 et  — = — ) 

S,  R, 


S1  V R 2 J [B  2 


u ± U-i 


Pour  une  surface  grise  Siéchangeant  du  rayonnement  avec  N surfaces  grises,  le  flux  radiatif  net  s’écrit: 

N 


= Y 

r,net,i  / 


^ r,net : i->j 


iogivi  cof  iiHi  a 


ijijij.  lesuise 
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_M°i-Ji_  y 

Vr,net,i  ± _ £ ^ 


Ji-Jj 


Si  X Ei 


7=1  Si  x Fij 


(33) 


L’échange  radiatif  entre  trois  surfaces  grises  S±,  S2  et  S3  constituant  une  enceinte  fermée  peut  être 
représenté  par  le  schéma  suivant  : 


1 £1  1 1 — £2 


Par  exemple,  le  flux  total  émis  par  5X  est  donné  par: 

3 


0 


= y 

,net,l  / 


0 


r,net:  1 -*/ 


7=1 


0, 


r,net,  1 


M°~Ji  _ V // 

1 - £1  ' ~ Z,  1 

5^  x ^ F ij 


> Exemples  d’application:  (voir  TD,  série  n°l ) 
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ANNEXE  1 


TABLE:  VALEURS  DE  LA  FONCTION  F0_àt(AT ),  FRACTION  DE  L ’ EMITTANCE  DU  CORPS 
NOIR  CORRESPONDANT  AUX  LONGUEURS  D’ONDE  INFERIEURES  A UNE  VALEUR  A (POUR 
UNE  TEMPERATURE  T DONNEE) 


AT 

(pin. K) 

0 

20 

40 

60 

80 

AT 

(pin. K) 

0 

20 

ou  200,,> 
ou  2000l+> 

40 

ou  4001'1 
ou  40007*' 

60 

ou  600' 
ou  6000l+> 

80 

ou  800,,> 
ou  8000l+> 

500 

0.0000 

0,0000 

0,0000 

0,0000 

0,0000 

6200 

0,7541 

0,7556 

0,7572 

0,7587 

0,7603 

600 

0.0000 

0,0000 

0,0000 

0,0000 

0,0000 

6300 

0,7618 

0,7633 

0,7648 

0,7662 

0,7677 

700 

0.0000 

0,0000 

0,0000 

0,0000 

0,0000 

6400 

0,7692 

0,7706 

0,7721 

0,7735 

0,7749 

800 

0.0000 

0,0000 

0,0000 

0,0000 

0,0000 

6500 

0,7763 

0,7777 

0,7791 

0,7804 

0,7818 

900 

0.0001 

0,0001 

0,0001 

0,0001 

0,0002 

6600 

0,7831 

0,7845 

0,7858 

0,7871 

0,7884 

1000 

0.0003 

0,0004 

0,0004 

0,0005 

0,0007 

6700 

0,7897 

0,7910 

0,7923 

0,7936 

0,7948 

1100 

0.0004 

0,0010 

0,0013 

0,0015 

0,0018 

6800 

0,7961 

0,7973 

0,7985 

0,7998 

0,8010 

1200 

0.0021 

0,0024 

0,0028 

0,0033 

0,0037 

6900 

0,8022 

0,8034 

0,8045 

0,8057 

0,8069 

1300 

0.0043 

0,0049 

0,0055 

0,0062 

0,0069 

7000 

0,8080 

0,8092 

0,8103 

0,8115 

0,8126 

1400 

0.0076 

0,0086 

0,0096 

0,0106 

0,0117 

7100 

0,8137 

0,8148 

0,8159 

0,8170 

0,8181 

1500 

0.0128 

0,0140 

0,0153 

0,0167 

0,0182 

7200 

0,8191 

0,8202 

0,8213 

0,8223 

0,8234 

1600 

0.0197 

0,0213 

0,0230 

0,0247 

0,0266 

7300 

0,8244 

0,8254 

0,8264 

0,8275 

0,8285 

1700 

0.0285 

0,0305 

0,0326 

0,0347 

0,0370 

7400 

0,8295 

0,8304 

0,8314 

0,8324 

0,8334 

1800 

0.0393 

0.0417 

0,0442 

0,0467 

0,0494 

7500 

0,8343 

0,8353 

0,8362 

0,8372 

0,8381 

1900 

0.0521 

0.0549 

0,0577 

0,0606 

0,0636 

7600 

0,8390 

0,8399 

0,8409 

0,8418 

0,8427 

2000 

0.0667 

0.0698 

0,0730 

0,0763 

0,0796 

7700 

0,8436 

0,8444 

0,8453 

0,8462 

0,8471 

2100 

0.0630 

0.0865 

0,0900 

0,0936 

0,0972 

7800 

0,8479 

0,8488 

0,8496 

0,8505 

0,8513 

2200 

0.1009 

0.1045 

0,1084 

0,1122 

0,1161 

7900 

0,8521 

0,8530 

0,8538 

0,8546 

0,8554 

2300 

0.1200 

0,1240 

0,1280 

0,1320 

0,1361 

8000 

0,8562 

0,8570 

0,8578 

0,8586 

0,8594 

2400 

0.1402 

0,1444 

0,1486 

0,1528 

0,1571 

8100 

0,8601 

0,8609 

0,8617 

0,8624 

0,8632 

2500 

0.1613 

0,1656 

0,1700 

0,1743 

0,1787 

8200 

0,8639 

0,8647 

0,8654 

0,8661 

0,8669 

2600 

0.1831 

0.1875 

0,1920 

0,1964 

0,2009 

8300 

0,8676 

0,8683 

0,8690 

0,8697 

0,8704 

2700 

0.2053 

0.2098 

0,2143 

0,2188 

0,2234 

8400 

0,8711 

0,8718 

0,8725 

0,8732 

0,8738 

2800 

0.2279 

0.2324 

0,2369 

0,2415 

0,2460 

8500 

0,8745 

0,8752 

0,8759 

0,8765 

0,8772 

2900 

0.2506 

0.2551 

0,2596 

0,2642 

0,2687 

8600 

0,8778 

0,8785 

0,8791 

0,8797 

0,8804 

3000 

0.2732 

0.2778 

0,2823 

0,2868 

0,2913 

8700 

0,8810 

0,8816 

0,8822 

0,8829 

0,8835 

3100 

0.2958 

0.3003 

0,3047 

0,3092 

0,3137 

8800 

0,8841 

0,8847 

0,8853 

0,8859 

0,8865 

3200 

0.3181 

0.3225 

0,3269 

0,3313 

0,3357 

8900 

0,8871 

0,8877 

0,8882 

0,8888 

0,8894 

3300 

0.3401 

0,3445 

0,3488 

0,3531 

0,3574 

9000 

0,8899 

0,8905 

0,8911 

0,8916 

0,8922 

3400 

0.3617 

0.3660 

0,3703 

0,3745 

0,3787 

9100 

0,8927 

0,8933 

0,8938 

0,8943 

0,8949 

3500 

0.3829 

0.3871 

0,3912 

0,3954 

0,3995 

9200 

0,8954 

0,8959 

0,8965 

0,8970 

0,8975 

3600 

0.4036 

0.4077 

0,4117 

0,4158 

0,4198 

9300 

0,8980 

0,8985 

0,8990 

0,8995 

0,9000 

3700 

0.4238 

0.4277 

0,4317 

0,4356 

0,4395 

9400 

0,9005 

0,9010 

0,9015 

0,9020 

0,9025 

3800 

0.4434 

0.4472 

0,4511 

0,4549 

0,4585 

9500 

0,9030 

0,9035 

0,9039 

0,9044 

0,9049 

3900 

0.4624 

0,4661 

0,4699 

0,4736 

0,4772 

9600 

0,9054 

0,9058 

0,9063 

0,9067 

0,9072 

4000 

0.4809 

0.4845 

0,4881 

0,4917 

0,4952 

9700 

0,9076 

0,9081 

0,9085 

0,9090 

0,9094 

4100 

0.4987 

0,5022 

0,5057 

0,5092 

0,5126 

9800 

0,9099 

0,9103 

0,9107 

0,9112 

0,9116 

4200 

0.5160 

0.5194 

0,5227 

0,5261 

0,5294 

9900 

0,9120 

0,9124 

0,9129 

0,9133 

0,9137 

4300 

0.5327 

0.5359 

0,5392 

0,5424 

0,5456 

10000 

0,9141 

0,9181 

0,9218 

0,9253 

0,9287 

4400 

0.5488 

0.5519 

0,5551 

0,5582 

0,5612 

11000 

0,9318 

0,9347 

0,9375 

0,9401 

0,9426 

4500 

0.5643 

0.5673 

0,5703 

0,5733 

0,5763 

12000 

0,9450 

0,9472 

0,9493 

0,9513 

0,9532 

4600 

0.5793 

0.5822 

0,5851 

0,5880 

0,5908 

13000 

0,9550 

0,9567 

0,9584 

0,9599 

0,9614 

4700 

0.5937 

0,5965 

0,5993 

0,6020 

0,6048 

14000 

0,9628 

0,9641 

0,9654 

0,9666 

0,9678 

4800 

0.6075 

0,6102 

0,6129 

0,6156 

0,6182 

15000 

0,9689 

0,9699 

0,9709 

0,9719 

0,9728 

4900 

0.6209 

0,6235 

0,6261 

0,6286 

0,6312 

16000 

0,9737 

0,9745 

0,9753 

0,9761 

0,9769 

5000 

0.6337 

0,6362 

0,6387 

0,6412 

0,6436 

17000 

0,9776 

0,9783 

0,9789 

0,9796 

0,9802 

5100 

0.6461 

0,6485 

0,6509 

0,6532 

0,6556 

18000 

0,9807 

0,9813 

0,9818 

0,9824 

0,9829 

5200 

0.6579 

0,6603 

0,6625 

0,6648 

0,6671 

19000 

0,9833 

0,9838 

0,9842 

0,9847 

0,9851 

5300 

0.6693 

0.6716 

0,6738 

0,6760 

0,6782 

20000 

0,9855 

0,9888 

0,9912 

0,9929 

0,9942 

5400 

0.6803 

0.6825 

0,6845 

0,6867 

0,6888 

30000 

0,9952 

0,9960 

0,9966 

0,9971 

0,9975 

5500 

0.6909 

0,6929 

0,6950 

0,6970 

0,6990 

40000 

0,9978 

0,9981 

0,9983 

0,9985 

0,9987 

5600 

0.7010 

0,7030 

0,7049 

0,7069 

0,7088 

50000 

0,9988 

0,9989 

0,9990 

0,9991 

0,9992 

5700 

0.7107 

0,7126 

0,7145 

0,7164 

0,7183 

60000 

0,9993 

0,9993 

0,9994 

0,9994 

0,9995 

5800 

0.7201 

0,7219 

0,7238 

0,7256 

0,7273 

70000 

0,9995 

0,9995 

0,9996 

0,9996 

0,9996 

5900 

0.7291 

0,7309 

0,7326 

0,7343 

0,7361 

80000 

0,9996 

0,9997 

0,9997 

0,9997 

0,9997 

6000 

0,7378 

0,7395 

0,7411 

0,7428 

0,7444 

90000 

0,9997 

0,9997 

0,9997 

0,9998 

0,9998 

6100 

0,7461 

0,7477 

0,7493 

0,7509 

0,7525 

100000 

0,9998 

0,9998 

0,9998 

0,9998 

0,9998 

(*)  pour  les  valeurs  de  AT  comprises  entre  10000  et  20000.  (+)  pour  les  valeurs  de  AT  supérieures  à 20000. 
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Chapitre  3:  Conduction  thermique 


1.  Introduction: 

On  rappelle  que  la  conduction  est  le  mode  de  transfert  de  chaleur  existant  dans  un  milieu  matériel  donné 
sans  qu’il  y ait  déplacement  macroscopique  de  matière.  C’est  ce  qui  se  passe,  en  particulier,  dans  les 
milieux  solides  homogènes  opaques  (parois  en  béton,  parois  métalliques...),  transparents  ou  semi- 
transparents  (parois  vitrées  comme,  par  exemple,  les  fenêtres...).  La  conduction  a lieu  également  dans  les 
fluides  immobiles,  c’est  le  cas,  par  exemple,  de  couche  d’air  immobile  confinée  entre  deux  couches  de 
matériaux  de  construction. 

La  conduction  ne  peut  se  produire  que  s’il  existe  un  écart  de  température  ( déséquilibre  thermique),  c’est-à- 
dire  que  si  le  gradient  de  température  n’est  pas  nul.  Dans  le  cas  contraire,  le  milieu  est  en  équilibre 
thermique  et  aucun  transfert  de  chaleur  ne  peut  avoir  lieu.  Pour  que  ce  gradient  de  température  existe,  il 
faut  une  action  externe  au  système  (sources  thermiques  externes)  pour  pouvoir  maintenir  des  conditions  de 
température  ou  de  flux  sur  les  frontières  du  système:  conditions  aux  limites  thermiques. 

Afin  de  simplifier  la  résolution  des  problèmes  de  transfert  de  chaleur  par  conduction,  l’étude  sera  faite  en 
régime  permanent,  pour  lequel  la  température  en  tout  point  du  milieu  est  indépendante  du  temps,  le 
gradient  de  température  (déséquilibre  thermique)  est  entretenu  par  les  sources  thermiques  extérieures. 

2.  Equation  générale  de  la  conduction: 

Pour  un  solide  homogène  isotrope  dont  les  propriétés  caractéristiques:  À,  Cp  et  p sont  constantes  et  pour 
lequel  p est  indépendante  de  la  température,  l’équation  régissant  le  transfert  de  chaleur  par  conduction  est 
donnée  par  (démonstration  voir  annexe  1 ): 

p(M,  t)  1 dT 

A T{M,  t)  + — = - X — (M,  t)  (1) 

a CL  ut 


Où  : 

■ p est  la  puissance  thermique  produite  (ou  consommée),  par  unité  de  volume  du  milieu,  par  la 
source  (ou  le  puits)  de  chaleur  interne.  Dans  le  cas  général,  p est  fonction  de  la  position  du  point 
M,  de  la  température  et  du  temps:  p(M,  T,  t).p>0  — > source',  p<0  — > puits . 

On  peut  citer,  comme  exemples,  chaleur  produite  (ou  consommée)  lors  d’une  réaction  chimique,  chaleur 
produite  par  effet  Joule  par  les  résistances  électriques  contenues  dans  un  four  électrique  ou  un  radiateur 
électrique... 


■ a — est  la  diffusivité  thermique  du  matériau.  Elle  caractérise  la  vitesse  à laquelle  la 

pXCp 

chaleur  se  propage  par  conduction  dans  un  matériau.  Plus  cette  diffusivité  est  grande,  plus  le 
matériau  s’échauffe  ou  se  refroidit  rapidement. 

On  rappelle  que  X est  la  conductivité  thermique  du  matériau  qui  traduit  son  aptitude  à conduire  la 
chaleur. 

La  valeur  de  X permet  de  quantifier  le  pouvoir  isolant  de  chaque  matériau.  Plus  elle  est  faible,  plus  le 
matériau  sera  isolant,  plus  elle  est  élevée,  plus  le  matériau  sera  conducteur. 

Le  produit  de  la  masse  volumique  p du  matériau  par  sa  chaleur  massique  ou  spécifique  Cp  représente 
la  capacité  thermique  de  ce  matériau.  Elle  exprime  sa  capacité  à absorber  et  à stocker  de  la  chaleur  en 
son  sein. 

Dans  le  tableau  1 , placé  en  fin  de  ce  chapitre,  on  donne  les  valeurs  de  X,  p et  Cp  de  quelques  corps 
usuels. 
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> Formes  particulières  de  l’équation  de  la  conduction 
On  rencontre  très  fréquemment  les  formes  particulières  suivantes  de  l’équation  de  la  conduction: 

milieu  avec  sources  (puits)  internes  de  chaleur,  en  régime  permanent. 

p 

AT  + — =0  (2)  Equation  de  Poisson 

A 

milieu  sans  sources  ni  puits  internes  de  chaleur,  en  régime  permanent. 

AT  = 0 (3)  Equation  de  Laplace 

milieu  sans  sources  ni  puits  internes  de  chaleur,  en  régime  variable : 

1 dT 

AT  = — x — - (4)  Equation  de  Fourier 

Ci  ot 

Comme  on  s’intéresse,  dans  ce  cours,  uniquement  au  transfert  de  chaleur  permanent  et  unidirectionnel 
dans  un  milieu  homogène  et  isotrope,  les  équations  à retenir  sont  celles  de  Laplace  et  de  Poisson : 


■ Equation  de  Laplace  : 

d2T 


(5)  ( Géométries  planes) 

(6)  ( Géométries  cylindriques) 


1 d 


= 0 


(7)  ( Géométries  sphériques) 


■ Equation  de  Poisson: 


d2T  p 
dxï  + J = ° 


(8)  ( Géométries  planes) 


1 9 1 

-x  — 

r dr 1 

t dT\ 

rx3“ 
v dr) 

+ 

II 

o 

(9) 

(Géométries  cylindriques) 

1 d 
r2  X dr 

( 2 97 

V X dr. 

+ 

II 

o 

(10) 

(Géométries  sphériques) 

3.  Condition  initiale  et  conditions  aux  limites: 

Pour  résoudre  l’équation  de  la  conduction,  il  faut  lui  associer  des  conditions  spécifiques  à chaque 
problème: 

la  condition  initiale  qui  donne  la  distribution  de  température  dans  tout  le  milieu  considéré,  à 
l’instant  initial,  t=0:  elle  utilisée  dans  le  cas  de  transfert  de  chaleur  en  régime  variable  ; 

les  conditions  aux  limites  qui  précisent  comment  la  chaleur  peut  traverser  les  frontières  du  milieu. 
Généralement,  on  distingue  trois  types  de  conditions  aux  limites  thermiques. 

3.1.  Condition  de  Direchlet  ou  de  lere  espèce: 

On  se  donne  en  chaque  point  de  la  surface  frontière  (D)  du  corps,  la  valeur  de  la  température  (température 
imposée): 

T(M)  = Cste  V M G (2)  (11) 

(25  est  donc  une  surface  isotherme. 


Exemple:  Cas  d’un  corps  en  contact  avec  un  thermostat. 


3.2.  Condition  de  Neumann  ou  de  2eme  espèce : 


dT 


La  densité  de  flux  de  chaleur,  c’est-à-dire  la  dérivée  normale  de  la  température  — (M)  est  connue  en  tout 


point  M de  la  surface  frontière  (. 
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-Ax  — (M)  = /(M)  V M E (?)  (12) 

on 

Exemple:  Corps  ayant  une  surface  adiabatique. 

3.3.  Condition  de  Fourier  ou  de  3eme  espèce  (condition  mixte): 

La  densité  de  flux  de  chaleur  sur  la  surface  frontière  (?)  est  contrôlée  par  un  coefficient  d’échange 
thermique  h connu  et  par  la  température  Tœ  du  fluide  environnant  loin  de  la  surface  frontière  du  système: 

dT 

-Ax—(M)  = hx(T(M)-Tœ)  V M e (2)  (13) 

on 

qui  s’écrit  encore: 

A dT 

T(M ) + — X — (M)  = T œ = Cste  V M E (?) 

h dn 

C’est-à-dire,  la  combinaison  linéaire  de  la  température  et  de  sa  dérivée  normale  est  connue  sur  la  surface 
frontière  du  corps. 

Exemple:  Cas  d’une  paroi  solide  en  contact  avec  un  fluide. 

3.4.  Contact  thermique: 

D’autres  conditions  aux  limites  peuvent  être  écrites  pour  décrire  le  passage  de  la  chaleur  d’un  corps  à un 
autre,  de  nature  différente,  à travers  leurs  surfaces  frontières  mises  en  contact  direct.  Même  si  les  surfaces 
sont  considérées  comme  lisses  la  chaleur  ne  transite  pas  nécessairement  avec  facilité  d’un  corps  à l’autre  et 
on  est  amené  à parler  de  contacts  thermiques  parfaits  ou  imparfaits. 


Soient  deux  corps  de  conductivités  thermiques  respectives  Ai  et  A2,  et  M un  point  appartenant  à leur 
surface  frontière  commune  (S).  Pour  chacun  des  corps,  les  flux  de  chaleur,  au  point  M,  s’écrivent: 

( dTx 

q>1  = —A\  x 


<p2  — —Ai  x 


dri! 

dT2 

dn7 


On  exprime  généralement  la  continuité  du  flux  thermique  au  contact  des  surfaces:  conservation  du 
flux  de  chaleur.  En  ce  qui  concerne  la  température,  il  existe  une  discontinuité,  un  saut,  lorsque  le 
contact  n’est  pas  parfait. 


3.4.1.  Contact  thermique  parfait: 


dT  1 dTz 

r,(M)  = t2(m> 


M E (S) 


(14) 


Fisure  1:  Contact  thermique  parfait 
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3.4.2.  Contact  thermique  imparfait: 

La  liaison  thermique  à l’interface  s’exprime: 

par  l’égalité  des  densités  de  flux  thermique  de  part  et  d’autre  de  l’interface: 

dT i dT2 

<p(M)  = -Ai  x (M)  = -A2x—± (M)  M g (5)  (15  - a) 

Oîlj  OTL2 


par  la  proportionnalité  du  saut  de  température  AT  = T 1 — T 2 au  flux  thermique: 

T1(M)-T2(M)  = rx  <p(M)  Me  (S)  (15-6) 


Le  coefficient  de  proportionnalité  r représente  la  résistance  thermique  de  contact  qui  sera  précisée  dans  la 
suite  de  ce  chapitre.  Elle  s’exprime  en  m2.°C/W  (ou  en  m2.°C.  h/Kcal)  et  dépend  de  la  rugosité  des 
surfaces,  de  leurs  propriétés  mécaniques,  de  la  pression  de  serrage  et  de  la  présence  de  gaz  interstitiel. 


<■> 

8 

Fisure  2:  Contact  thermique  imparfait 


4.  Analogie  électrique,  notion  de  résistance  thermique: 

En  régime  permanent,  le  flux  de  chaleur  traversant  un  mur,  sans  sources  internes  de  chaleur,  d’épaisseur  e, 
de  conductivité  A,  de  grandes  dimensions  transversales  (L  et  l)  et  dont  les  faces  extrêmes  sont  maintenues  à 
des  températures  uniformes  et  constantes  J;  et  T2.  ( T i>T2),  est  donné  par  la  relation: 


0= Ax S x 


T , - T,  T T7 


(16) 


AxS 


Où  est  S est  l’aire  de  la  surface  traversée  par  la  chaleur:  S=Lxl. 

L et  / sont  choisies  grandes  devant  e pour  satisfaire  à l’hypothèse  du  transfert  unidirectionnel. 


La  deuxième  forme  de  la  relation  (16)  est  analogue  à la  loi  d’Ohm  en  électricité  qui  définit  l’intensité  du 
courant  comme  le  rapport  de  la  différence  de  potentiel  électrique  à la  résistance  électrique.  La  température 


apparaît  ainsi  comme  un  potentiel  thermique  et  le  terme  — apparaît  comme  la  résistance  thermique  d’un 

AS 


mur  simple  d’épaisseur  e,  de  conductivité  A et  de  surface  S. 


Red  ~ 


AxS 


(17) 


Red 

V www-  - 

O 


Ti 


0=i 


Ti-T2 
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Rcd  est  la  résistance  au  transfert  de  chaleur  par  conduction  entre  deux  surfaces  isothermes.  Elle 
s’exprime  en  °C/W  ou  en  °C.h/Kcal  et  caractérise  la  difficulté  pour  réaliser  un  transfert  de  chaleur.  Pour  un 
même  écart  de  température  entre  deux  surfaces,  la  puissance  thermique  transférée  sera  d’autant  plus 
importante  que  la  résistance  thermique  Rcd  sera  faible. 

Si  l’on  multiplie  Rcd  par  S on  obtient  la  résistance  thermique  par  unité  de  surface  notée  rcd: 

rcd  = | (18) 

Elle  s’exprime  en  m2.°C/W  ou  en  m2.°C.h/Kcal. 


Si  l’une  des  parois  du  mur  (par  exemple  celle  portée  à T t)  est  en  contact  avec  un  fluide  à la  température  Tm 
le  flux  thermique  échangé  par  convection  s’écrit: 

0cv  = hcxSx  (T1  -Tœ)=  Tl~Tœ  = Tl~  Tœ  (19) 

"ci; 

hc  x S 

Donc  : 


Rcv  hr  X S 


(20) 


Rcv  est  la  résistance  thermique  relative  à l’échange  de  chaleur  par  convection  entre  la  paroi  du  mur  et  le 
fluide  environnant. 


Si  la  paroi  et  le  fluide  échange  de  la  chaleur  par  rayonnement,  le  flux  radiatif  net  échangé  entre  eux  s’écrit: 


0 , . ...  = 5 , x J2"  x air j —Tt) 

net / paroi— > fluide  1 12  \ 1 00  / 


(21) 


Cette  expression  peut  être  linéarisée  en  fonction  de  Ti-Tx  dans  le  cas  de  faibles  écarts  de  température 

(jj)  —7^1  <10(Pc). 

,net\paroi-> fluide  ~ tlr  X S X (T i — (22) 

avec: 


hr  -47rn  a x 


T +T 

1 J T1œ 


(23-a) 


Par  analogie  avec  la  convection,  le  paramètre  hr  ainsi  défini  est  appelé  coefficient  d’échange  thermique 
par  rayonnement. 

Dans  le  cas  contraire  >10(PC ),  le  coefficient  h,  se  calcule  en  utilisant  la  relation  suivante: 

Fuxo(Tf-T*)  , 2 2s  , s 

K = +TJ  (23-b) 


La  relation  (22)  peut  être  écrite  comme  suit: 

*t*r,net-.paroi-+ fluide  tlr  X S X (T ^ T m) 


T i -Tn 


hr  x S 


Donc  : 


Rr  K x S 


(25) 


Tx-Tœ 

Rr 


(24) 


Rr  est  la  résistance  thermique  relative  à l’échange  de  chaleur  par  rayonnement  entre  la  paroi  du  mur  et 
le  fluide  environnant. 

Pour  déterminer  une  résistance  thermique  totale,  comme  en  électricité,  on  sait  que: 
en  série,  deux  ou  plusieurs  résistances  s’ajoutent, 
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en  parallèle,  les  inverses  des  résistances  thenniques  s’ajoutent. 


Rj  R ->  Rj  Rn  Rthea 

= Mi \àr^- 

0 0 


Figure  4: 

n 

Rth,eq  = ^ ' Ri  (26) 


f=l 


Mpw 


T, 


0 


Ri 


0 


*WH Un 

Mpw 


Figure  5: 


-=t- 

Rth,eq  Ri 


(27) 


5.  Applications: 

On  se  place  en  régime  permanent  et  on  suppose  que  les  milieux  conductifs  sont  dépourvus  de  sources  (ou 
de  puits)  internes  et  que  la  chaleur  s’y  propage  suivant  une  seule  direction. 

5.1.  Mur  simple,  mur  composite: 

5.1.1.  Mur  simple: 

Soit  un  mur  simple  d’épaisseur  e de  conductivité  A et  dont  les  deux  faces  planes  sont  portées  aux 
températures  J)  et  T2  ( T î>T2 ).  On  se  propose  de  déterminer  la  distribution  de  température  dans  ce  mur  ainsi 
que  le  flux  thermique  le  traversant. 


On  affaire  à un  mur  simple  à surfaces  frontières  isothermes,  c’est-à-dire  soumises  à des  conditions  aux 
limites  de  Direchlet  ou  de  Ie"  espèce.  Le  système  fondamental  donnant  la  distribution  de  température  à 
l’intérieur  du  mur  s’écrit: 


( d?T 


T(x  = 0)  = T1 
<T(x  = e)  = T2 


(28) 


En  intégrant  ce  système,  on  obtient: 

T i - T 2 

T(x)  = xx  + Tx  pour  0<x<e  (29) 

e 

C’est  la  distribution  de  température  à l’intérieur  du  mur. 


La  densité  de  flux  thermique  traversant  ce  mur  s’écrit,  d’après  la  loi  de  Fourier. 

dT  T i - T 2 

<p(x)  = -Ax  — = Ax— (30) 
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Le  flux  thermique  traversant  une  surface  d’aire  S du  mur  s’écrit: 

j1  y 

0 = S x <p(x ) = Ax  — - 


Soit,  en  introduisant  la  résistance  thermique  de  conduction: 

T i - T 2 

0=-^ (32) 


R 


cd 


Avec  : 


Rcd  ~ 


AxS 


(31) 


5.1.2.  Mur  composite  formé  de  murs  plans  homogènes,  en  série: 

Il  est  facile  de  généraliser  les  résultats  obtenus  pour  un  mur  simple  à un  mur  composite  constitué  de 
plusieurs  couches  de  matériaux  accolées  d’épaisseur  e,  et  de  conductivité  Ai  en  contact  thermique  parfait. 

Si  les  surfaces  extrêmes  sont  maintenues  aux  températures  constantes  7/  et  1' le  flux  de  chaleur 
traversant  le  mur  composite  est  donné  par: 


Tj  ~ l’w+l 
Zf=l  Rcdi 


T\  ~ Tn+ 1 

Rcd,eq 


(33) 


Si  les  surfaces  extrêmes  sont  en  contact  avec  deux  fluides  de  températures  Txi  et  {figure  6),  le  flux  de 
chaleur  traversant  le  mur  composite  est  donné  par: 


0 = 


Rcv  1 


T ool  _ T 002 
+ Xf=l  Rcdi  + Rcv 2 


(34) 


avec: 


Rcdi  ~ — — : résistance  thermique  de  conduction  du  matériau  i (e,-,  Aj)  ; 


Rcv  1 
Rcv  2 


A,xS 

1 

hixS 

1 

tl2XS 


: résistance  thermique  de  convection  du  côté  du  fluide  (7); 
: résistance  thermique  de  convection  du  côté  du  fluide  (2). 
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5.1.3.  Mur  composite  constitué  de  plusieurs  murs  plans  homogènes,  en  parallèle: 

Supposons  maintenant  que  différents  éléments  solides  soient  juxtaposés  par  bandes  et  que  la 
température  soit  uniforme  sur  chacune  de  leurs  faces  (figure  7).  L’écart  de  température  TrT2  est  donc 
le  même  pour  chacun  des  éléments  traversé  respectivement  par  les  flux  thermiques  ( £>;,...  0h ...  et  0,. 


Si  Rj,  R2,. . .,  Ri,...  et  Rn  représentent  les  résistances  thermiques  de  conduction  de  chacun  des  murs 
simples,  alors  les  flux  les  traversant  respectivement  02,  02,.. . 0,...  et  0,  sont  donnés  par: 


0,= 


T -T 

± 1 L 2 

R, 


02  = 


T -T  T -T 

1 2 . /T,  _ 1 1 l2 


R, 


R: 


et  0„ 


Ti 


Le  flux  thermique  total  traversant  le  mur  composite: 


0>™=2>,=(r,-r2£;— = 


b i Ti  ~T2 


i=l 


R, 


theq 


La  résistance  thermique  équivalente  s’écrit: 


D _ 

theq  n y 


x,;  z 

i=l 


f AS> 


i=l  \ c J 


i-1 


(35) 


D'où:  R[heq  = — (36) 

Zfcs,) 

i=l 

Remarque: 

Si  les  différents  murs  simples  en  parallèle  n’ont  pas  la  même  épaisseur,  le  raisonnement  précédent 
s’applique  à condition  de  négliger  les  échanges  thermiques  par  les  faces  latérales  des  murs  juxtaposés. 


5.2.  Conduite  cylindrique,  paroi  cylindrique  composite: 

5.2.1.  Conduite  cylindrique: 

Soit  une  conduite  cylindrique  de  très  grande  longueur  L par  rapport  aux  rayons  intérieur  r2  et  extérieur  r2.  et 
dont  les  faces  interne  et  externe  sont  soumises  à des  températures  uniformes  et  constantes  Tj  et  T2  (T/>T2). 
Le  matériau  constituant  la  conduite  est  de  conductivité  thermique  A. 


Les  isothermes  sont  des  surfaces  cylindriques  coaxiales,  le  gradient  de  température  est  radial  (la 
température  n’est  fonction  que  de  r).  Le  système  fondamental  s’écrit  donc: 

1 d ( dT\ 

X^(rX^)  = ° 


T(r  = rx)  = T1 
Tir  = r2)  = T2 


(37) 


et  dont  la  résolution  aboutit  à l’expression  suivante  de  T(r): 


WWW. 
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T(r ) = 7\  + 


Ln 


(R) 


x Ln  ^ ^ pour  r-±  < r < r2  (38) 


On  déduit  de  la  loi  de  Fourier  la  densité  de  flux,  <p: 


T1  - T 2 1 

(p(r)  = A X t zr-r  X 


in 


m 


(39) 


Le  flux  thermique  traversant  chaque  surface  cylindrique  de  rayon  r et  de  longueur  L est: 

0=  S(r ) x (p{r ) = (2.7 z.r.L)  x cp(r ) 


D’où: 


0 = 2.n.  L.  A x 


Ln 


(R) 


(40) 


Enfin,  on  déduit  l’expression  de  la  résistance  thermique  de  conduction  d’une  conduite  cylindrique'. 


Ln 


Rcd  - 


(R) 


2 x 7T  x L x A 


(41) 


5.2.2.  Paroi  cylindrique  composite  (conduites  cylindriques  coaxiales  accolées): 

Le  flux  thermique  traversant  la  paroi  des  conduites  cylindriques  coaxiales  accolées  de  conductivités 
thermiques  différentes  se  rencontre  dans  de  nombreuses  installations  industrielles.  Un  tuyau  isolé 
transportant  un  fluide  chaud  et  exposé  à une  ambiance  froide  est  l’un  des  problèmes  types  (figure  8).  Si  le 
tuyau  est  relativement  long,  la  direction  de  la  chaleur  s’écoulant  à travers  sa  paroi  serait  radiale.  On 
supposera  un  contact  thermique  parfait  entre  les  conduites. 


Fluide  froid  Tx2,  h2 

Figure  8:  conduites  cylindriques  accolées 
La  conduite  (l)=(Ai,rh  r2)  constitue  une  première  résistance  au  transfert  thermique  Rc,ii- 

(r_F 
UV 


^cdl  ~ 


2 x n x L x A1 

La  conduite  (2)=(A2,r2,  r2)  constitue  une  deuxième  résistance  au  transfert  thermique  Rc,p. 

fîV 


Rcdi  — 


2xnx  Lx  A2 

La  résistance  thermique  de  conduction  équivalente  est: 
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D’où: 


(42) 


Le  flux  thermique  traversant  la  paroi  cylindrique  composite  est: 

T i - T3 

0 = -4 1 (43) 


fl 


cd,eq 


Si  et  T3  sont  connues, 
ou  bien: 


0 = 


CO  1 * CO  2 


CO  1 * CO  2 


Rev  1 Rcd,eq  Rev 2 


fl 


th,eq 


Si  Lœl  et  rœ2  sont  connues, 
avec: 


Rcvl 


Il  1 / ,V  1 X (2.7T.  Ti-L) 


(44) 


fl,.„i  = = -:  résistance  thermique  de  convection  entre  le  fluide  chaud  et  la  paroi  interne 

CV1  h!XS ! h1x(2.n.r1.L)  M 1 

de  la  conduite  (7); 

1 1 

Rri, 7 = = résistance  thermique  de  convection  entre  la  paroi  externe  de  la  conduite 

cv L h2xS3  h2x(2.n.r3.L)  M 1 

(2)  et  le  fluide  froid. 


Ce  résultat  peut  être  généralisé  pour  une  paroi  cylindrique  composite  constituée  de  n conduites  coaxiales: 

T ,TJ\  — T oo2 


0 = 


R 


th,eq 


avec: 


R th,eq  Rcvl  Rcd,eq  “I"  Rev 2 Rev  1 


" Ln 


+ 


L 

i= 1 


ft?) 


2xn  x Lx  Ai 


+ R 


cv  2 


(45) 


5.3.  Sphère  creuse,  paroi  sphérique  composite : 

5.3.1.  Sphère  creuse: 

Considérons  une  sphère  creuse  de  rayons  intérieur  rt  et  extérieur  r2,  et  dont  les  surfaces  frontières  sont 
soumises  à des  températures  uniformes  et  constantes  7)  et  T2.  Les  isothermes  sont  des  surfaces  sphériques 
concentriques.  Le  transfert  thermique  est  donc  radial.  Le  système  fondamental  s’écrit: 


1 d 


x 


dT 

dr 


H 


T(r  = rf)  = 7\ 
Tir  = r2)  = T2 


(46) 


L’intégration  de  ce  système  conduit  à la  solution  suivante: 

1 _ _1 

r(r)  = 7\  + {T 2 - Tf)  x (n  < r < r2)  (47) 


On  déduit  de  la  loi  de  Fourier  la  densité  de  flux: 

T!-T2 


<p(r)  = 2 x 


WWW. 


x 


(48) 
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et  ensuite  le  flux  thermique  0 traversant  chaque  surface  sphérique  de  rayon  r: 

0=  5(r)  x (p{r ) = (4. 7T.  r2)  x <p(r) 


D’où: 


0 = 4.71.  A x 


Lri  r 2 J 


(49) 


Enfin,  on  déduit  l’expression  de  la  résistance  thermique  de  conduction  d’une  sphère  creuse: 

1 1 


r 1 r2  r2  — rx 

4.7 1.  A 4.n.A.r1.r2 


(50) 


5.3.2.  Paroi  sphérique  composite  (sphères  concentriques  accolées): 

Le  flux  thermique  traversant  la  paroi  de  deux  sphères  concentriques  accolées  (figure  8)  S (r1,r2,Ai)  et  S 
(r2,r3,A2)  est  donné  par: 

Ti~T3 

0=^ 1 (49) 

Kcd,eq 

Si  et  T 3 sont  connues, 
ou  bien: 


0 = 


T ooi  - T , 


00  2 


T 00 1 - T, 


00  2 


R 


cvl 


+ Rcd,eq  + R 


cvï 


R 


th,eq 


(51) 


Si  T ooi  et  Tœ2  sont  connues, 
où: 


R , — R ,,  + R 

cd.eq  cal  cal 

composite; 


4n 


V2  "''Z 
y^irir2 


+ - 


^2r2r3  J 


est  la  résistance  thermique  de  conduction 


de  la  paroi 


Rcvl  — = r js  est  la  résistance  thermique  de  convection  entre  la  paroi  interne  de  la  sphère  Si 

h j S j hjX  r 4 jtr  1 I 

et  un  fluide  environnant; 


R , = = y est  la  résistance  thermique  de  convection  entre  la  paroi  externe  de  la  sphère 

^2^3  h 2 A2  ^3  ) 

S et  un  fluide  environnant. 


www . ne su 


Figure  9:  Sphères  concentriques  accolé \ 
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5.4.  Epaisseur  critique  d’un  revêtement  isolant: 

Si  dans  le  cas  d’une  paroi  plane  la  moindre  épaisseur  d’isolant  entraîne  une  réduction  des  transferts 
thermiques,  il  n’en  est  pas  de  même  pour  une  conduite  cylindrique  ou  une  paroi  sphérique.  En  fait, 
l’addition  d’une  couche  isolante  à la  surface  extérieure  d’une  conduite  ne  réduit  pas  toujours  les  pertes 
thermiques  vers  le  milieu  ambiant  mais  contribue  parfois  à l’augmentation  de  ces  pertes,  jusqu’à  atteindre 
un  maximum  de  pertes  pour  une  épaisseur  donnée  appelée  épaisseur  critique  d’isolation. 

5.4.1.  Cas  d’une  conduite  cylindrique: 

On  se  propose  d’isoler  une  conduite  cylindrique,  de  conductivité  thermique  X,  de  diamètres  interne  Dj  et 
externe  D2  et  de  longueur  L ( L»D1  et  D2),  par  un  manchon  isolant  d’épaisseur  x et  de  conductivité  A,-.  On 
suppose  que  la  surface  interne  de  la  conduite  est  en  contact  avec  un  fluide  chaud  à la  température  Txj, 
coefficient  d’échange  thermique  h],  et  que  la  surface  externe  du  manchon  est  en  contact  avec  un  fluide 
froid  à la  température  Tx2,  coefficient  d’échange  thermique  h2.  Quelle  est  l’influence  de  l’épaisseur  du 
manchon  isolant  sur  les  pertes  thermiques? 


Figure  10:  Conduite  cylindrique  isolée. 


Les  pertes  thermiques  à travers  la  paroi  de  la  conduite  nue  sont  données  par: 

T . -T 


0 - 

v o 

=>  <Z>  = 


00  1 00 2 


R J + R , + R -y 

cvl  cd  cv2 


T -T 

00 1 *052 


(53) 


1 1 _ 

1 Lu 

rd)1Lh1  2n2L 


rD  A 


\Di  j 


+ 


7tD2Lh2 


Les  pertes  thermiques  à travers  la  paroi  de  la  conduite  isolée  sont  données  par: 


0 = 


T -T 

* 1 1 oo2 


R .+/?,+/?,.  R , 

cvl  cd  edi  cv2 


0 = 


T -T 

col  002 


1 1 

H x Ln 


rD  ^ 


7d)1Lh1  2xAL 


\Di  j 


+ 


Ou  encore: 


0{x)  = 


2 nk:  L 


x Ijî 


1+2. 


D 


2 J 


n{p , +2x)lJi2 


2tiLxRi,  h xLn 

tri 


1 + 2. 


(54) 


+ 


■ 't  l D2  J (D2+2*L 
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avec  : R.,  — — 1 — — xLn 

th  7d)1Ui1  2 7i AL 


\Di  J 


Calculons  la  dérivée  de  0 par  rapport  à l’épaisseur  x du  manchon: 


(10 

dx 


~\2 


2Æx/?„  H xLu 

th  A. 


1 + 2, 


D 


2 J 


(d2  +2x)i2 


X{D2+2xYlX 


f x ' 

D?  +2x—2.— 

2 ]hj 


(55) 


On  voit  que  cette  dérivée  est  susceptible  d’être  positive  dans  certaines  conditions. 

A. 

En  effet,  si  Z)  >2 — la  dérivée  de  0 est  toujours  négative  quel  que  soit  x et  l’augmentation  de 

k2 

l’épaisseur  de  l’isolant  réduit  les  pertes  thermiques. 

A.  /I,  I), 

Par  contre,  si  Z),  <2—  il  existe  une  valeur  xc  = — , épaisseur  critique  de  l’isolant,  pour  laquelle 


h2  2 


les  pertes  thermiques  sont  maximales.  Ceci  est  dû  au  fait  que  l’augmentation  de  la  surface  d’échange 
externe  (diminution  de  la  résistance  de  convection)  provoquée  par  l’addition  de  la  couche  isolante 
l’emporte  sur  l’augmentation  de  la  résistance  de  conduction.  Lorsque  x varie  de  0 à xc,  les  pertes 
augmentent  de  0„=0(x=O)  à 0max=0(x=xc),  puis  décroissent  en  passant  par  0O  obtenu  pour  une  épaisseur 
Xo- (figure  10).  Pour  réduire  les  pertes  dans  ce  cas,  il  faudrait  dépasser  Xo  tel  que: 


2tiLxR„  + — xLii 

th  A. 


' xn ' 

1+2. — - 

V D2J 


+ 


Soit,  en  remplaçant  0()  par  son  expression: 

Lu 


f x ^ 

1 + 2 — — 

V D2  J 


A. 

-2  — x 

h , 


D 2 +2X0  J 


{d2+2x0)i2 

(56) 
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(W) 


Figure  11:  Variation  des  pertes  thermiques  traversant  une  conduite,  de  petit  diamètre 
externe  ( D2=12,5  mm),  en  fonction  de  l’épaisseur  de  l’isolant. 


H 


h. 


(W) 
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g\ 

n\ 

n\ 

i\  \ 


Isolant  de  conductivité  Àn=h2x(D2/2)=l,5  W/m.°C 
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\\  \ Magnésie  (fc=0, 078  W/m.°C) 

_l\  \ V Laine  de  verre  (A=0,045  W/i?i.°C) 
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\ 


. \ Polyuréthane  (A=0, 035  W/m.°C) 

W V v 


Isolant  de  conductivité  A<h,x(D,/2) 

- _ 0" 


0.05 

Epaisseur  de  l’isolant  x 


Figure  12:  Variation  des  pertes  thermiques  traversant  une  conduite,  de  grand  diamètre 
J externe  (Di^dO  cn^%pfonëtiém  de  l’épaisseur  de  l’isolant. 
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5.4.2.  Cas  d’une  paroi  sphérique: 

À.  D 

Un  raisonnement  analogue  conduit  à x =2— , épaisseur  critique  de  la  couche  sphérique  isolante. 

h 2 2 

Par  conséquent,  pour  réaliser  une  isolation  thermique  efficace  il  faut  choisir  un  isolant  de  conductivité  A* 

D 2 

telle  que:  Â.  < Ag,  avecA0  — h2  . 


8.  Ailettes,  surfaces  ailetées: 

Dans  de  nombreuses  situations  pratiques  on  est  amené  à augmenter  le  transfert  de  chaleur  entre  une  surface 
donnée  (plane,  cylindrique  ou  sphérique)  et  un  fluide  environnant  en  fixant  des  ailettes  sur  cette  surface. 

Les  ailettes  sont  utilisées  lorsqu’il  faut  échanger  une  puissance  thermique  importante  dans  un 
encombrement  réduit.  Comme  exemples,  on  peut  citer  le  refroidissement  des  composants  électroniques,  les 
radiateurs  d’automobile,  les  échangeurs  industriels... 

D’une  manière  générale,  l’emploi  des  ailettes  est  avantageux  pour  les  gaz  qui  sont  caractérisés  par  un 
coefficient  d’échange  convectif  h faible,  et  moins  utile  pour  les  liquides  pour  lesquels  h est  grand.  Les 
ailettes  minces,  étroites  et  rapprochées  sont  meilleures,  au  point  de  vue  transfert  de  chaleur,  que  les  ailettes 
épaisses  et  espacées.  Cependant,  on  est  limité  par  les  pertes  de  charge  qui  augmentent  si  l’on  diminue  trop 
l’espacement  des  ailettes.  Une  ailette  est  d’autant  plus  performante  que  sa  conductivité  thermique  est 
élevée. 

Le  choix  d’une  forme  géométrique  convenable  d’une  ailette  exige  un  compromis  entre  le  prix 
(raisonnable),  l’encombrement  (réduit),  pertes  de  charge  (faibles)  et  le  transfert  de  chaleur  (important). 


6.  Conduction  en  régime  permanent  avec  sources  internes  de  chaleur: 

Les  systèmes  avec  sources  (ou  puits)  de  chaleur  sont  rencontrés  dans  de  nombreuses  applications 
techniques.  Les  exemples  types  sont  les  réacteurs  chimiques  et  nucléaires,  les  chaudières,  les  radiateurs 
électriques... 

Dans  ce  chapitre,  on  ne  considérera  que  quelques  cas  simples  pour  lesquels  la  production  interne  de  la 
chaleur  est  uniforme  (plaque  plane,  cylindre,  sphère).  L’équation  générale  de  la  conduction  se  réduit  à celle 
de  Poisson  : 

P 

Ar  + ^ = o 

A 

6.1.  Plaque  plane  avec  sources  internes  uniformes: 

Considérons  une  plaque  plane  et  épaisse  dont  les  surfaces  frontières  sont  portées  à des  températures 
uniformes  et  constantes,  et  dans  laquelle  la  chaleur  est  générée  uniformément  (figure  13).  La  longueur  et  la 
largeur  de  la  plaque  sont  supposées  grandes  devant  son  épaisseur  pour  que  le  transfert  de  chaleur  y soit 
unidirectionnel.  En  régime  permanent,  l’équation  de  la  chaleur  s’écrit: 
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dont  l’intégration  conduit  à la  solution  suivante: 


T(r)=-^r2  +C1Ln(r)+C2 


où  Ci  et  C2  seront  déterminées  à l’aide  des  conditions  aux  limites  associées  au  problème. 

Comme  il  s’agit  d’un  cylindre  plein  de  rayon  R dont  la  surface  latérale  est  soumise  à une  température 
uniforme  et  constante  Ta.,  les  conditions  aux  limites  s’écrivent  comme  suit: 


dr 


= 0 


pas  d' échappement  de  chaleur  au  centre  du  cylindre 


r=0 


[T(r  = R)  = T0 


Cj=0 


c2  =T0  +-^-R2 
2 ° 41 


Par  conséquent,  la  distribution  de  température  dans  le  cylindre  est  donnée  par: 


~r2YTo  ( 0<r<R ) 


(62) 


La  température  maximale  est  localisée  sur  l’axe  du  cylindre  et  égale  à: 


Tnwx  =T(rm  =0)=T0  + 


pRz 


4A 


(63) 


6.3.  Sphère  pleine  avec  sources  internes  uniformes: 

Soit  une  sphère  pleine  de  rayon  R,  de  conductivité  A,  siège  d’une  production  uniforme  de  chaleur  p, 
plongée  dans  un  fluide  de  température  Tx  avec  un  coefficient  d’échange  thermique  h.  Le  système  présente 
une  symétrie  de  révolution,  le  transfert  thermique  est  donc  radial.  Le  système  fondamental,  composé  de 
l’équation  de  conduction  et  des  conditions  aux  limites,  s’écrit: 


1 d 


r2  dr 


2dT > 


V 


dr 


-A— 

dr 


= 0 (pas  d' échappement  de  chaleur  au  centre  de  la  sphère) 


r=0 


dr 


= h{T{r  = R)-T„ ) 


r=R 


puissance  produite  à l'intérieur  de  la  sphère 
et  libérée  au  fluide  par  deux  mécanismes 
agissant  en  parallèle  : convection  + rayonnement 


L’intégration  de  ce  système  donne  la  solution  suivante: 


T(r)  = To0  + 


pR- 


r 


6A 


r p 
1-—  +—R 

2 3h 


Rz 


(O  <r  <R)  (64) 


De  la  loi  de  Fourier  on  déduit  le  flux  thermique: 


&(r)=-7tr3  p=V(r)p 


La  température  sur  la  surface  de  la  sphère  s’obtient  en  écrivant: 

0(r)=vp=Iis(t(r)-tx) 


-n R3p  = hS(T(R)-T„)  (s=4Æ2) 


www,  n 


(65) 


. com 
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Tableau  1:  Propriétés  Thermoühxsiaues  de  certains  corps 


Corps 

Point  de 
fusion 

Température 

Densité 

Chaleur 

spécifique 

Conductivité 

thermique 

Notation  & Unité 

Tf(°C) 

T(°C) 

P (Kg/m3) 

Cp  (J/Kg.  °C) 

A (W/m.°C) 

Métaux,  alliages  et  céramiques 

Aluminium  pur 

660 

20 

2700 

900 

204 

Alumine  ( AUO 3) 


Argent  pur 


Cuivre  pur 


Laiton  (70%  Cu,  30%  Zn) 


Fer  pur 


Acier  inox  (15%  Cr,  10%  Ni) 


Acier  inox  (18%  Cr,  8%  Ni) 


Plomb 


Tungstène 


Alumine  (AUO 3) 


Magnésie 


20 


20 


400 


20 


400 


20 


20 


20 


20 


20 


20 


38 


Matériaux  isolants 


10525 


8954 


249 


29 


407 


386 


363 


111 


147 


73 


20 


16,3 


35 


176 


29 


0,067 


Laine  de  roche 


Laine  de  verre 


Laine  minérale 


Liège  expansé 


Polyuréthane  (mousse  dure) 


PVC  (mousse  rigide) 


Polystyrène  expansé 


20 

55 

880 

0,038 

20 

135 

880 

0,041 

20 

10 

880 

0,045 

20 

40 

880 

6,035 

20 

15 

839 

0,04 

20 

120 

2100 

0,044 

20 

32 

1300 

0,03 

20 

50 

1360 

0,035 

20 

85 

1300 

0,045 

20 

30 

1300 

0,031 

20 

40 

1300 

0,041 

20 

28 

1300 

0,037 

Corps 

Température 

Densité 

Chaleur  spécifique 

Conductivité  thermique 

Notation  & Unité 

T(°C) 

p (Kg/m3) 

Cp  (J/Kg.  °C) 

A (W/m.  °C) 

Liquides 

Eau 

20 

1000 

4176,6 

0,6 

Huile  légère 


Glycérine 


Mercure 


(Tf=-38,9°C  et  Te=357°C) 


www . iesui 


962 


883 


865 


910 


865 


1270 


1250 


1547 


1397 


1216 


13707 


13409 


! 13168  | 135,6  | 12,63  I 

s etudiant . com 


4185 


1674 


1841,4 


1841,4 


2134 


2318,5 


2511 


875 


934,5 


1022 


141,5 


136,5 


0,675 


0,16 


0,15 


0,132 


0,128 


0,286 


0,283 


0,0675 


0,0727 


0,0675 


9,76 


10,38 
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Sodium 

100 

927 

1384 

85,7 

(Tf=97,6°C  et  Te=882°C) 

500 

833 

1262 

66,3 

Potassium 

100 

824 

813 

48  ! 

(Tf=63,2°C  et  Te=757°C) 

500 

721 

762 

37,5 

NaK  (56%  K,  44%  Na) 

100 

880 

1064 

22,8  j 

(Tf=5,7°C  et  Te=812°C) 

500 

881 

982 

27,1 

Lithium 

204 

509 

4365 

46,37 

(Tf=179°C  et  Te=1317°C) 

538 

476 

4171 

30,45 

Bismuth 

315 

10011 

144 

16,4 

(Tf=271°C  et  Te=1477°C ) 

538 

9739 

155 

15,6  \ 

Plomb 

371 

10540 

159 

16,1 

(Tf=327°C  et  Te=1737°C) 

704 

10140 

155 

14,9 

Gaz 

-128,9 

2,462 

1001 

0,0137 

Air 

38 

1,136 

1005 

0,023 

538 

0,438 

1101 

0,0574  j 

100 

0,595 

1885 

0,025  \ 

Vapeur  d’eau 

538 

0,271 

2134 

0,067  ! 

1093 

0,16 

2511 

0,131 

-128,9 

2,719 

911 

0,0137 

Oxygène 

93 

1,07 

931 

0,0311 

538 

0,482 

1058 

0,0628 

-240 

0,724 

10300 

0,0242 

Hydrogène 

93 

0,0673 

14396 

0,21  ! 

538 

0,0304 

14689 

0,36 

-128,9 

0,338 

5230 

0,09 

Hélium 

93 

0,134 

5230 

0,17 

538 

0,0604 

5230 

0,315  1 

-17,8 

2,116 

795 

0,0133  \ 

Gaz  carbonique 

93 

1,472 

913 

0,022 

538 

0,668 

1169 

0,0549 
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ANNEXE  1 

EQUATION  DE  BASE  DE  LA  CONDUCTION 

Considérons  un  solide  opaque,  isotrope,  de  volume  ( VJ,  de  chaleur  spécifique  cp.  de  conductivité  thermique 
A,  et  dont  la  température  T varie  dans  l’espace  et  dans  le  temps,  T(M,t).  Avec,  en  plus,  des  sources  (puits) 
internes  de  chaleur  générant  (consommant)  une  puissance  thermique  par  unité  de  volume p(M,t)  (W/ni). 

1.  Forme  générale  de  l’équation  de  la  conduction: 
y Bilan  énergétique: 

On  applique  le  principe  de  conservation  de  l’énergie  à un  volume  élémentaire  dV  de  surface  frontière  dS 
contenu  ( V),  pendant  l’intervalle  du  temps  dt  (figure  1). 


Le  bilan  énergétique  peut  s’écrire,  littéralement,  sous  la  forme: 

Chaleur  qui  entre  dans  (dV)  pendant  dt-Chaleur  qui  sort  pendant  dt+Chaleur  dégagée  ( consommée ) 
dans  (dV)  par  les  sources  (ou  puits)  internes  pendant  dt+Chaleur  due  à la  variation  de  l’énergie  interne 

pendant  dt. 

Ou  encore: 

Entrée  — Sortie  + Génération  (ou  Dégradation ) = Accumulation 


Sous  une  forme  algébrique: 


(<PxdSx  + <PydSy  + <PzdSz  )~  (<Px+dxdS  x + 9. y+dyd^  y + <Pz+dzdS z )+  p(M ,t)dV  = fjC p ^~dV 


Entrée 


Sortie 


Génération  ( dégradation  ) 


Accumulation 


( i ) 


(<Px  <Px+dx  )dS  dx  \ ?y  P }• +‘iy  )ds  dy  + <Pz+dz\lS,dz  + p(M,t)dV  = pC  —dV 
dx  x dy  y dz  p dt 


S(px  , 8<Py  , d(pz 


dx  dy  dz 


I dV  + p(M,t)dV  = pC„  — dV 
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Or:  tp  = 


<Py 


, vecteur  densité  de  flux  de  chaleur. 


L’équation  du  bilan  devient: 


(-  div(p+ p(M,t)}lV  = pCp  —dV 

dt 

d‘V  étant  un  volume  arbitraire  de  ( V),  on  a donc  en  tout  point  de  ( V)\ 

dT 

-divtp+ p(M,t)  = pC  p — 

dt 

Compte  tenu  de  la  loi  de  Fourier  (tp  = —AgradT  ),  cette  équation  devient: 

div{Agrad(T(M ,t))j+  p(M,t)  = pC  p ^ (M,t ) 
C’est  l’équation  de  la  chaleur  dans  un  solide  isotrope. 


(2) 


2.  Formes  particulières  de  l’équation  de  la  conduction: 

Pour  un  solide  homogène  dont  les  propriétés  caractéristiques:  A,  C/>  et  p sont  constantes  et  pour  lequel  p est 
indépendante  de  la  température,  on  obtient  l’équation  linéaire  classique  de  la  chaleur: 

Ar+jKMP=Ler  (3) 

A a dt 

où  a = est  la  diffusivité  thermique  du  matériau. 

pC  p 


3.  Expressions  analytiques  de  l’équation  de  la  chaleur: 

En  introduisant  les  expressions  du  Laplacien  A en  coordonnées  cartésiennes  (x,y,z),  cylindriques  (r,9,z)  et 
sphériques  (r,6,tp),  on  obtient  respectivement: 


l_d_ 
r2  8r 


82T  d2T  ô2T  p 


- + 


1 dT 


dx2  Ôy2  dz2 
l_d_(  dT)  1 d2T  d2T 
rôry  dr  J r2  dO2  dz2 


Ut  A J Cl 

+ -[x/y/z/t)=-— 

A a ôt 


(4) 


A a dt 


(5) 


ÔT\ 


- +- 


r 


dr  J rswQ  d6\ 


swQ — + 


ô2t 


P(  O A ldT 

I 2 ■ 2 o a 2 5 v / 0,(p,t) 

dd  J r sw  9 dtp  A a dt 


(6) 
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Chapitre  4:  Convection  Thermique 


1.  Généralités  et  définitions: 

1.1.  Différentes  formes  de  convection: 

La  convection  est  un  transport  de  chaleur  dû  à des  mouvements  macroscopiques  du  fluide.  Elle 
intervient  en  particulier  dans  les  échanges  thermiques  entre  une  paroi  solide  et  un  fluide. 

On  distingue  deux  formes  de  convection:  la  convection  libre  et  la  convection  forcée . 

La  distinction  se  fait  sur  l’origine  du  mouvement  permettant  le  transfert  d’énergie  d’une  région  à une 
autre. 

Lorsque  ce  mouvement  est  provoqué  par  les  seules  différences  de  températures,  donc  de  densité,  on  dit 
que  la  convection  est  libre  ou  naturelle. 

Lorsque  le  fluide  est  mis  en  mouvement  par  un  dispositif  mécanique  (ventilateur,  pompe...),  on  dit 
que  la  convection  efforcée. 

Lorsque  les  deux  formes  de  convection  coexistent,  on  dit  que  la  convection  est  mixte. 

1.2.  Régimes  d’écoulement: 

L’étude  de  la  convection  est  intimement  liée  à celle  des  régimes  d’écoulement  de  fluide.  Ainsi,  dans 
l’écoulement  d’un  fluide  le  long  d’une  paroi  solide  deux  régimes  principaux  peuvent  être  observés:  ou 
bien  le  fluide  s’écoule  par  filets  parallèles  à la  paroi,  le  régime  est  alors  laminaire,  ou  bien  à cet 
écoulement  se  trouve  superposé  une  agitation  correspondant  à des  vitesses  normales  à la  paroi  qui,  en 
un  point  donné,  sont  constamment  variables  dans  le  temps  mais  restent  petites  devant  les  vitesses 
générales  de  l’écoulement  (valeurs  moyennes  des  composantes  de  la  vitesse  pendant  un  temps  très 
long).  On  dit  alors  que  le  régime  est  turbulent. 

11  existe  un  troisième  régime  d’écoulement  instable  qui  apparaît  entre  les  régimes  laminaire  et 
turbulent  et  auquel  on  donne  le  nom  de  régime  intermédiaire  ou  troublé. 


Régime 

Régime  laminaire  troublé 

k >k >k 

• • é 

\ 

v Repos 

Figure  1: 


Régime  turbulent 


* y 


moy 

► 


Re 


1.3.  Couches  limites: 


On  considère  une  paroi  solide  plane  (figure  2)  qu’on  suppose  normale  au  plan  de  la  figure  soit  OX  sa 
trace  sur  le  plan;  elle  est  balayée  par  un  courant  de  fluide  parallèle  à l’axe  OX  dont  la  vitesse  à l’infini 
est  constante,  soit  Uœ.  Le  bord  d’attaque  normale  en  0 se  projette  en  O.  A cause  de  la  viscosité,  la 
vitesse  est  nulle  sur  la  paroi  et  elle  croît  régulièrement  de  O à Uœ.  quand  on  se  déplace  dans  le  sens  de 
OY.  En  réalité,  les  variations  de  la  vitesse  sont  pratiquement  localisées  dans  une  mince  couche  fluide 
au  voisinage  de  la  plaque  dite  couche  limite.  Le  calcul  attribue  à cette  couche  une  section  parabolique 


(figure  2);  son  épaisseur  à l’abscisse  x étant 


Sr 


, v est  la  viscosité  cinématique,  v 


est  un  coefficient  purement  numérique. 


1 
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Il  se  forme  également  depuis  le  bord  d’attaque  de  la  paroi  une  couche  limite  thermique  dans  laquelle 
sont  localisées  les  variations  de  température  (figure  3). 


Loin  du  bord  d’attaque  de  la  paroi  solide,  les  deux  couches  se  superposent.  On  obtient  ainsi  une 
couche  fluide  laminaire  dont  l’épaisseur  dépend  des  propriétés  physiques  du  fluide  et  de  la  nature  de 
l’écoulement.  Au  sein  de  cette  couche  de  fluide,  on  admet  qu’il  n’y  a aucun  mélange  de  matière  et  que 
la  chaleur  ne  peut  s’y  transmettre  que  par  conduction. 

La  conductivité  thermique  des  fluides  étant  faibles  par  rapport  à celles  des  solides,  la  couche 
laminaire,  bien  que  mince,  oppose  une  résistance  importante  au  passage  de  la  chaleur.  Il  en  résulte 
dans  cette  région  une  variation  rapide  de  la  température. 

En  dehors  de  la  couche  laminaire,  au  sein  du  fluide,  la  chaleur  se  transmet  par  mélange  des  particules 
de  fluide,  provoquant  ainsi  une  égalisation  rapide  de  la  température:  on  parle  ainsi  de  la  température 
du  mélange  du  fluide,  T «, 


Fluide  circulant  à U 0r=C'“’  et  T 0r=Cs" 


C c r q 

Cj  n r»  CO 


wz. 


C C r 

ne  Cj 


Couche  limite 
laminaire  mince 
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Dans  le  cas  d’un  écoulement  de  fluide  à l’intérieur  d’une  conduite  cylindrique,  les  deux  couches 
limites  thermique  et  dynamique  s’élargissent  à mesure  qu’on  s’éloigne  du  bord  d’attaque  (entrée  de  la 
conduite)  jusqu’à  occuper  tout  l’intérieur  de  la  conduite,  on  dit  alors  que  les  régimes  dynamique  et 
thermique  sont  établis  (figure  5). 


Fluide  U, 


Figure  5: 


Zone  où  les  régimes  dynamique 
et  thermique  sont  établis 

mmmMMmmMMmmmmmMmmmm. 

Zone  où  les  régimes  dynamique 
et  thermique  ne  sont  pas  établis  U=0 


2.  Calcul  du  coefficient  d'échange  thermique  par  convection,  h : 

Considérons  un  fluide  à la  température  du  mélange  Tœ,  en  contact  avec  une  paroi  solide  d’aire  S et  de 
température  Tp  (Tp  < Tœ)  (figure  3).  Le  flux  de  conduction  à travers  la  couche  laminaire  peut  être 
exprimé  par  la  loi  de  Fourier  : 

a) 

e 

où  e est  l’épaisseur  de  la  couche  laminaire  et  X la  conductivité  thermique  du  fluide. 

Cependant,  cette  relation  n’est  pas  utilisable  en  pratique,  car  on  ne  connaît  généralement  pas 
l’épaisseur  e de  la  couche  laminaire. 

X 

On  exprime  alors  le  flux  par  la  loi  de  Newton,  en  posant  h = — : 

e 

0 = hxSx(Tœ-Tp)  (2) 
et  on  cherche  à évaluer  le  coefficient  h. 


Pour  définir  le  coefficient  de  convection,  on  distinguera  deux  cas:  celui  de  la  convection  à l’intérieur 
d’une  conduite  cylindrique  et  celui  d’une  surface  plongée  dans  un  fluide  indéfini.  On  va  introduire  les 
notions  essentielles  concernant  la  convection  forcée,  en  régime  permanent,  en  analysant  un  cas 
particulier  correspondant  à l’écoulement  de  fluide  dans  une  conduite  cylindrique.  Ces  notions  sont 
valables  pour  les  autres  géométries.  Le  cas  de  la  convection  naturelle  sera  également  traité  à la  fin  de 
ce  chapitre. 


2.1.  Convection  forcée  à l’intérieur  d’une  conduite  cylindrique'. 

Le  calcul  du  coefficient  de  convection  h n’est  pas  aisé.  L’expérience  montre  qu’en  l’absence  de 
convection  naturelle,  h n’est  pas  influencé  par  la  différence  de  température  (Tœ  —Tp),  mais  qu’il 
dépend  des  sept  grandeurs  suivantes: 

Vm:  vitesse  moyenne  du  fluide, 
p:  masse  volumique  du  fluide, 
cp : chaleur  massique  ou  spécifique  du  fluide, 
ju:  viscosité  dynamique  du  fluide, 

X:  conductivité  thermique  du  fluide, 

D y.  diamètre  intérieur  de  la  conduite, 

x : abscisse  de  la  section  considérée,  par  rapport  à l’entrée  de  la  conduite. 
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On  peut  donc  écrire: 

h = f(ym,  p,  Cp,  P,  A,  Dh  x)  (3) 

Il  est  pratiquement  impossible  de  déterminer  expérimentalement  l’influence  de  chacune  des  grandeurs 
sur  h,  c’est-à-dire,  de  connaître  la  fonction  /.  11  faudrait,  pour  cela,  ne  faire  varier  qu’une  seule 
grandeur  en  maintenant  les  autres  constantes.  Or,  cela  n’est  pas  possible  pour  un  nombre  aussi  élevé 
de  variables,  en  raison  d’inévitables  interactions. 

Cette  complexité  des  phénomènes  de  convection  rend  nécessaire  l’utilisation  de  techniques  générales 
permettant  de  limiter  le  nombre  de  paramètres  d’un  problème.  Parmi  ces  techniques  on  peut  citer 
/ ’ analyse  dimensionnelle. 

2.1.1.  Analyse  dimensionnelle: 

Cette  méthode  stipule  que  la  relation  / liant  entre  elles  huit  variables  (h,  Vm,  p,  cp,  p,  À,,  D,,  x)  qui 
s’expriment  au  moyen  de  quatre  unités  fondamentales  ( M : masse  (Kg),  L : longueur  (m),  t : temps  (s), 
T : température  ( °C ))  peut  être  réécrite  sous  la  forme  d’une  autre  fonction  F qui  lie  8 — 

Nbre  de  variables 

4 =4  nombres  adimensionnels. 

Nbre  d'unités 

Unités  fondamentales  : 


M : masse  (Kg) 

L : longueur  (m) 

t : temps  (s) 

T : température  (°C) 

Dimensions  des  variables: 


[h]  :M  L2  f3.  L2.  T'=M  f3  T1  (W/m2.  °C) 

[V,„]  : LJ1  (m/s) 

[p]  ML3  (Kg/m3) 

[Cp\  : L2  f2  T1  (J/Kg.  °C) 

[p\  : M L1  f1  (Kg/m.s)  ou  (Pa.s) 

[A]  : M L f3  T1  (W/m.  °C) 

[Di]  : L (m) 

[x]  • L (m) 


La  relation  f(h,  Vm,  p,  cp,  p,  A,  ü,,  x)  = 0 peut  être  remplacée  par: 


hx  Di 


V „ (PXVmX  Di\ak  (P  X CV  ( X \ 

= 2>x( — ^ — ) x(— ) XU) 

k=  1 


X \ck 


(4) 


Cette  relation  est  appelée  corrélation  entre  nombres  sans  dimension.  Pour  préciser  sa  forme,  il  faut 
déterminer  les  constantes  Ak  et  les  exposants  (ak,  bk,  ck)  en  menant  une  étude  expérimentale.  Les 
nombres  sans  dimension  ainsi  créés  (4)  sont  classiques. 

hD 

Nu  — — — : nombre  de  Nusselt,  il  caractérise  l’échange  thermique  entre  le  fluide  et  la  paroi. 


WV5  ^ 


Re  = 


A 

pVmDi 


P 


nombre  de  Reynolds,  il  caractérise  le  régime  d’écoulement  du  fluide 


(circulation  forcée): 

si  Re  < 2000,  le  régime  d’écoulement  est  laminaire, 

si  2000  < Re  < 3000,  le  régime  d’écoulement  est  intermédiaire  ou  troublé, 

si  Re  > 3000,  le  régime  d’écoulement  est  turbulent. 

^ q gn  i ^ in 
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> 


> 


Pr  — 


p 

"77' 


nombre  de  Prandtl,  il  caractérise  les  propriétés  dynamiques  et  thermiques  du 


fluide.  Il  mesure  le  rapport  de  la  diffusivité  des  quantités  de  mouvement  à la  diffusivité 
thermique  du  fluide. 

-7  ; abscisse  réduite. 

Di 


En  résumé,  les  travaux  expérimentaux  étudiant  le  transfert  de  chaleur  par  convection  forcée  dans  une 
situation  donnée  (géométrie,  fluide...)  fournissent  leurs  résultats  sous  forme  de  corrélations 
mathématiques: 

n ^ 

Nu  = F (fle,  Pr, ^-)  = ^ Ak  x x x (5) 

1 fc= î 1 


Les  nombres  Re,  Pr  et  sont  calculables  dès  que  sont  connus  les  paramètres  du  problème.  On  peut 

alors  en  déduire  h en  calculant  d’abord  le  nombre  de  Nusselt,  Nu,  par  la  formule  (5),  puis  h = ~~~ 
(voir  TD). 


N. B.:  pour  calculer  la  valeur  numérique  des  nombres  adimensionnels,  il  faut  bien  veiller  à,  toujours, 
utiliser  un  système  cohérent  d’unités,  par  exemple,  le  système  international. 

2.1.2.  Principales  corrélations  en  convection  forcée: 

2.I.2.I.  Régime  d’écoulement  laminaire  (Re<Rec=2000): 

> Pour  des  conduites  très  longues,  le  régime  thermique  est  établi: 

dans  le  cas  où  la  température  de  la  paroi  de  la  conduite  est  constante  si  — > 0,035. Re  .Pr  , on 
utilisera: 

Nu  = 3,  66  (6) 

dans  le  cas  où  le  flux  traversant  la  paroi  de  la  conduite  est  constant  si  ^>0,l.Re  .Pr , on 
utilisera: 

Nu  = 4,  3 6 (7) 


Ces  valeurs  limites  de  Nu,  obtenues  analytiquement,  dépendent  de  la  géométrie  de  la  section  de  la 
conduite. 


r Pour  des  conduites  de  longueur  finie  ( ^ < 0,035.  Re  . Pr  ou  ^ < 0,1.  Re  . Pr  ),  le  régime 
thermique  n’est  pas  établi,  on  pourra  utiliser: 


la  corrélation  de  Hausen  valable  pour  — < 0,05. Re  .Pr  : 


F 


0,0668. 


Nu  = 3,66  + 


Re .Pr . 


D 


V 


1 + 0,04. 


Re  .Pr . 


D 


(8)  ( température  constante  de  la  paroi  de  la  conduite ) 


ou 
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f 


0,0668. 


Nu  = 3,66  + 


Re  .Pr . 


D 


V 


f \0’14 


1 + 0,04. 


Re  .Pr 


2 

D \3 


P P J 


(9) 


P, 


— — : facteur  de  correction  pour  la  viscosité. 

MP 

Les  indices  met  p indiquent  que  la  viscosité  dynamique  du  fluide  est  mesurée  à la  température 
moyenne  du  fluide  (T fm)  respectivement  à la  température  de  la  paroi  solide  (Tp). 


la  corrélation  de  Sieder  et  Tate  valable  pour 

.0,14 


D 


Re.Pr.—  \>10: 
L 


Nu  =1,86. 
avec:  0,48  <Pr<  16700 


f D) 

i 

( \ 

U 

1 3 

i tu 

l LJ 

(10)  (température  constante  de  la  paroi  de  la  conduite ) 


2.I.2.2.  Régime  d’écoulement  turbulent  (Re >Re c=3 000 ) : 

> Pour  des  conduites  très  longues  ( — >60  ),  on  pourra  utiliser: 

la  corrélation  de  Colburn  pour  Pr  > 0,5  et  104<Re<12  X 104; 

Nu  = 0, 023  x Re0,8  x Pr033  (11) 


la  corrélation  de  Dittus-Boelter,  valable  pour  0, 7 < Pr  < 120  et  Re  > 104: 
Nu  = 0, 023  x Re0  8 x Pr11  (12) 


avec;  n = 0,  3 si  Tp  < T fm  (refroidissement  du  fluide) 
n = 0, 4 si  Tp  > T fm  (échauffement  du  fluide) 

(T p température  de  la  paroi  de  la  conduite,  T f m température  moyenne  du  fluide). 


la  corrélation  de  Séban-Shimazaki  pour  les  métaux  liquides  (Pr  <<  1),  si  la  température  de 
la  paroi  de  la  conduite  est  constante: 

Nu  = 4,8  + 0,025  x (fie  x Pr)0' ^ (13) 


la  corrélation  de  Lyon-Martinelli  valable  pour  les  métaux  liquides  à faible  nombre  de  Prandtl 
(Pr  « 1)  et  pour  Re  x Pr  > 50: 

Nu  = 7 + 0, 025  x (Re  x Pr)0'8  (14) 


Soulignons  que  les  corrélations  de  Colburn  et  de  Dittus-Boelter  (relations  (11)  et  (12))  doivent  être 
appliquées  pour  des  différences  de  température  \T p — T fm\  faibles  voire  modérées,  avec  toutes  les 
propriétés  du  fluide  évaluées  à T fm. 


Pour  les  grands  écarts  de  température  Tp  — T fm,  il  est  recommandé  d’utiliser  la  corrélation  suivante 
dite  de  Sieder  et  Tate,  valable  pour  Re  > 104  et  0,  7 < Pr  < 16700: 


Nu  = 0, 027  x fle0'8  x Pr^  x ( — 


0,14 


(15) 


r Pour  des  conduites  de  longueur  finie  ( — < 60),  on  pourra  utiliser  la  corrélation: 


D 


WV5  ^ 


îûQHi  cof  iiHi  a 
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Nu  = Nu(Colburn ) x 


(16) 


> Pour  une  conduite  en  serpentin  (A  pas  du  serpentin  et  Dt  son  diamètre  interne),  on  pourra 
utiliser  la  corrélation: 


Nu  = Nu(Colburn ) x 


1 + 1,35  x 


D; 


(17) 


Les  grandeurs  physiques  À,  jt i,  p,  etc...  dont  dépendent  les  nombres  adimensionnels  sont  prises 

T fe  + ^ fs 

habituellement  à la  température  moyenne  du  fluide  Tfin  — , entre  l’entrée  ( x = 0)  et  la  sortie 

(x  = L ) de  la  conduite.  Cependant,  si  la  paroi  est  à une  température  Tp  très  différente  de  T fm  , il 


conviendra  d’introduire  une  autre  moyenne  Tm  = 


fin 


[ P 

— qui  est  la  température  du  film  laminaire. 


Pour  les  conduites  à sections  non  circulaires,  on  peut  utiliser  toutes  les  corrélations , à condition 
toutefois  de  remplacer  le  diamètre  interne  de  la  conduite  figurant  dans  l ’ expression  de  divers  nombres 
adimensionnels  par  le  diamètre  hydraulique  D h équivalen  t donné  par: 

4x5 

Dh  = (18) 

où  5 est  l’aire  de  la  section  de  la  veine  fluide  et  P le  périmètre  mouillé  interne. 


Exemples: 

> Conduite  à section  rectangulaire: 


Dh  = 


4 x ax  b 2 x axb 


2 x (a  + b)  a + b 


(19) 


Dh  = 


> Espace  annulaire: 
4 x (52  - 5X) 


P^  +P? 


— D-i  — D i 


(20) 


Figure  6: 


Figure  7 


> Espace  entre  deux  plans  parallèles: 

C’est  l’extension  de  la  formule  relative  à la  conduite  à section  rectangulaire  si  a -*  oo. 


Dh  = 2 x b (21) 


Figure  8 


b 


Fluide 


mmm  iesuisé-bndf^nt.  . roin 
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2.2.  Convection  forcée  à F extérieur  d’obstacles  solides 
2.2.1.  Cas  de  V écoulement  d’un  fluide  le  long  d’une  plaque  : 
a)  Ecoulement  laminaire  pour  Re  < 3 X 105  : 

> Pour  Pr  > 0,5,  on  utilisera  la  corrélation  suivante  : 

1 1 

Nu  = 0, 664  x Re 2 x Pr3  (22) 


> Pour  Pr  « 1 (métal  liquide,  on  utilisera  la  corrélation  suivante  : 


Avec  : 


Nu  = 1, 1 x yl (1  — Pr) 3 x (Re  x Pr)  (23) 

p x Vœ  x L 


Re  = 


P 


(24) 


Nu  = 


hx  L 

À 


(25) 


L étant  la  longueur  de  la  plaque  et  Vœ  la  vitesse  du  fluide  loin  de  la  paroi. 

b)  Ecoulement  turbulent  pour  Re  > 5 X 105  : 

> Pour  Pr  > 0,5,  on  utilisera  la  corrélation  suivante  : 

4 1 

Nu  = 0,035  x Pe5  x Pr3  (26) 

> Pour  Pr  « 1 (métal  liquide),  on  utilisera  la  corrélation  suivante  : 

Nu  = 0,  59  x (Re  x pr)0'61  (27) 


Pour  les  deux  régimes  d’écoulement,  les  propriétés  physiques  du  fluide  sont  prises  à la  température  du 
film  laminaire  Tm. 


2.2.2.  Cas  de  la  circulation  d’un  fluide  perpendiculairement  à l’axe  de  la  conduite: 

Pour  une  circulation  du  fluide  perpendiculaire  à la  conduite,  il  se  forme  un  sillage  en  aval  de  celle-ci 
dont  la  forme  et  la  dimension  dépendent  du  nombre  de  Reynolds , Re. 


Pour  calculer  le  nombre  de  Nusselt  Nu  on  utilisera: 

la  corrélation  de  Hilpert  si  le  fluide  est  un  gaz: 

Nu  = A x Rem 

avec: 

p x Vm  x De 
Re  = 

P 


(28) 

(29) 
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De  est  le  diamètre  extérieur  de  la  conduite  et  Vm  est  la  vitesse  moyenne  du  fluide  calculée  en  amont  de 
la  conduite  (loin  de  la  paroi). 

la  corrélation  de  Knudsen  et  Katz  si  le  fluide  est  un  liquide: 

Nu  = 1, 11  x A x Rem  x Pr0  31  (30) 

avec: 

Le  coefficient  A et  l’exposant  m étant  les  mêmes  pour  les  deux  corrélations;  leurs  valeurs  à utiliser 
sont  données  dans  le  tableau  suivant: 


Re 

A 

m 

Nu  {Gaz) 

1-4 

0,891 

0,330 

0,891-1,41 

4-40 

0,821 

0,385 

1,40-3,40 

40-4000 

0,615 

0,466 

3,43-29,34 

4000-40000 

0,174 

0,618 

29,3-121,5 

40000-400000 

0,024 

0,805 

121,6-776 

Pour  les  conduites  à section  non  circulaire,  on  utilisera  la  corrélation: 

1 

Nu  = A x Rem  x Pr3  (31) 

Les  valeurs  de  A et  m sont  données  dans  le  tableau  suivant  pour  différents  types  de  section  de  la 
conduite: 


Re 

A 

m 

V ^ 

- 1 

Q 

◄ ► 

5xl03-105 

0,102 

0,675 

M 

V a 

ly  De 

* t 

5x103-105 

0,246 

0,588 

H 

w 

< ► 

«5 

5x1 01- 1,95x1 04 

1,95  xlO4 -105 

0,160 

0,0385 

0,638 

0,782 

M 

Plaque 

U 

W De 

b 

verticale 

4000-1, 5xl04 

0,228 

0,731 

2.2.3.  Cas  de  la  circulation  d’un  fluide  perpendiculairement  à un  faisceau  de  conduites  cylindriques 
parallèles: 

L’écoulement  perpendiculaire  à des  rangées  de  tubes  a une  importance  considérable  tout 
particulièrement  dans  les  réchauffeurs  d’air  et  les  échangeurs  de  chaleur  tubulaires  ( voir  chapitre  3 du 
cours  des  machines  thermiques ). 
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Figure  10 -a:  Faisceau  aligné. 


Figure  1 0-b  : Faisceau  quinconcé. 


L’expérience  montre  que  la  disposition  de  tubes  en  quinconce  est  le  siège  d’une  grande  turbulence, 
donc  d’un  coefficient  de  convection  plus  élevé,  que  la  disposition  en  ligne. 


Pour  déterminer  la  valeur  du  nombre  de  Nusselt , Nu,  Colburn  a recommandé  la  corrélation  suivante: 

Nu  = C x fie06  x Pr0'33  (32) 


avec: 


| C —0,26  p oui'  un  faisceau  aligné 
IC  =0,33  pour  un  faisceau  quinconcé 


valable  pour  2 X 103  < Re  < 4 X 104. 


fie  = 


P * Vmax  x De 

P 


(33) 


Vmax  représente  la  vitesse  dans  la  section  la  plus  étroite. 


2.3.  Calcul  du  coefficient  d’échange  thermique  h en  convection  naturelle: 

Dans  la  convection  naturelle,  le  mouvement  de  fluide  est  dû  aux  variations  de  sa  masse  volumique 
provenant  des  échanges  de  chaleur  entre  le  fluide  et  la  paroi. 

Ces  échanges  jouent  un  rôle  en  pratique,  associés  souvent  à des  échanges  par  rayonnement.  Citons  par 
exemple: 

le  chauffage  domestique, 

les  pertes  par  les  parois  dans  les  installations  industrielles. 

2.3.1.  Description  du  phénomène: 

Considérons  une  plaque  plane  verticale  chaude  (dont  la  température  de  la  paroi  est  TP)  au  contact  d’un 
fluide  froid  (dont  la  température  loin  de  la  paroi  est  Tx).  Au  voisinage  de  la  plaque  existe  une  zone 
perturbée  appelée  couche  limite  de  convection  naturelle. 

2.3.2.  Résultats  de  l’analyse  dimensionnelle: 

Dans  le  cas  d’un  transfert  de  chaleur  par  convection  naturelle  le  long  d’une  plaque  plane,  le  coefficient 
de  convection  h dépend  des  caractéristiques  du  fluide:  X p,  p.  cP.  [3.  g,  de  la  paroi  caractérisée  par  la 
longueur  L,  et  de  l’écart  de  température  AT  aux  bornes  du  film  ce  que  l’on  peut  traduire  par  une 
relation  du  type: 

h = f{X,  p,  fj.,C  P,  f3,  g,  L,  AT  ) (34) 

Des  considérations  d’analyse  dimensionnelle,  semblables  à celles  exposées  dans  le  cas  de  la 
convection  forcée,  amènent  à définir  un  groupe  adimensionnel  appelé  nombre  de  Grashof. 

WWW . ’î&SUÎSgtn<ii  ant  rorn 
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ou: 


Gr  = 


L3.p2.g.jB.AT 

M2 


(35) 


L (m):  est  une  dimension  linéaire  permettant  de  calculer  la  surface  d’échange  (hauteur  de  la 
plaque  verticale  dans  l’exemple  précédent,  côté  d’un  carré,  largeur  d’un  rectangle,  diamètre 
d’un  cylindre,  etc.), 

P CK1)  : est  le  coefficient  de  dilatation  volumique  du  fluide  à pression  constante  (pour  l’air  et 


les  gaz:  f}  = 


T' 


aVeC  T',noy  = ~ X 


Tp  + 


f T +T  Y) 

co(  entrée  ) oo(  sortie  ) 

2 yy 


p (Kg/m. s)\  est  la  viscosité  dynamique  du  fluide, 
p (Kg/m3)  : est  la  masse  volumique  du  fluide, 
g:  est  l’accélération  de  la  pesanteur  ( g=9,81  m/s2), 

£T=lTrTj  (°C  ou  °K):  est  la  différence  de  température  entre  la  paroi  solide  et  le  fluide. 


Le  nombre  de  Grashof  Gr,  joue  le  même  rôle  en  convection  naturelle  que  le  nombre  de  Reynolds,  Re, 
en  convection  forcée.  Sa  valeur  conditionne  la  nature  laminaire  ou  turbulente  de  l’écoulement  dans  la 
zone  perturbée. 


Indiquons  que  le  nombre  sans  dimension: 

Ra  = Gr  x Pr  (36) 

s’appelle  nombre  de  Rayleigh  et  est  souvent  utilisé  dans  la  littérature. 

Si  l’on  reprend  l’exemple  de  la  plaque  verticale,  plus  la  hauteur  L de  celle-ci  augmente,  plus  la  valeur 
du  nombre  de  Grashof  croît. 

L’expérience  montre  que,  pour  les  gaz,  lorsque  ce  nombre  dépasse  109,  les  filets,  jusque  là  laminaires, 
présentent  des  ondulations  qui  s’amplifient  progressivement  et  deviennent  turbulents  dans  la  partie 
supérieure. 


Figure  11: 


En  convection  naturelle,  le  nombre  de  Nusselt,  Nu,  s’exprime  en  fonction  des  nombres  de  Grashof, 

Gr,  et  de  Prandtl,  Pr : 

Nu  — F(  Gr,Pr) 

Les  corrélations  sont,  en  général,  de  la  forme: 

Nu  = Cx  (Gr  x Pr)n  (37) 

www  - i^suisg'hiKii  ant*  corn 
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Les  grandeurs  physiques  sont  déterminées  à la  température  du  film  laminaire: 
f F t . t 'N’X 


T'  = — x 

moy  2 


Tp  + 


T + T 

oo(  entrée  ) co(  sortie  ) 


jj 


Les  valeurs  du  coefficient  C et  de  l’exposant  n sont  données  dans  le  tableau  suivant: 


Géométrie  et 
orientation  de  la  paroi 

Dimension 

caractéristique  pour  le 
calcul  de  Gr  et  de  Nu 

Coefficient  C 

Régime  laminaire 

n=l/4 

Régime  turbulent 

n=l/3 

Plaque  horizontale  chauffant 
vers  le  haut 

Largeur  1 

0,54 

( 1 05<(  Gr.Pr)<2  xl  O7) 

0,14 

(2x107<(  Gr.Pr)<4xl010) 

Plaque  horizontale  chauffant 
vers  le  bas 

Largeur  I 

0,27 

(3x1 05<(  Gr.  Pr)<3xl  O10) 

0,07 

(3xl010<(  Gr.  Pr)<1013) 

Plaque  verticale  ou  cylindre 
vertical  de  grand  diamètre 
par  rapport  à l ’ épaisseur  de 
la  zone  perturbée 

Hauteur  H ( ou  D ) 

0,59 

(1 04<(Gr.Pr)<l  O9) 

0,13 

(1 09<(Gr.Pr)<l  O13) 

Cylindre  horizontal 

Diamètre  extérieur  D 

0,53 

(1 03<(Gr.Pr)<l  O9) 

0,10 

(1 09<(Gr.Pr)<l  O13) 

Remarque'. 

La  valeur  de  l’exposant  n dépend  du  régime  de  convection. 

La  valeur  du  coefficient  C dépend  du  régime  de  convection  ainsi  que  de  la  géométrie  et  de 
l’orientation  de  la  paroi. 

Si  une  plaque  plus  chaude  que  le  fluide  environnant  a sa  face  chauffante  tournée  vers  le  haut,  le  fluide 
chaud,  plus  léger,  s’élève,  et  il  est  remplacé  facilement  par  du  fluide  froid;  si  au  contraire  la  face 
chauffante  est  tournée  vers  le  bas  le  fluide  chaud  n’est  que  difficilement  remplacé  par  du  fluide  froid; 
les  échanges  par  convection  devront  donc  être  plus  faibles  si  la  face  chauffante  est  tournée  vers  le  bas. 
C’est  bien  ce  que  montrent  les  valeurs  du  coefficient  C du  tableau  ci-dessus.  Ces  considérations 
doivent  être  inversées  si  la  plaque  est  plus  froide  que  le  fluide. 

On  utilisera  ainsi  la  même  corrélation  pour  une  face  chaude  vers  le  haut  que  pour  une  face  froide 
vers  le  bas. 


> Exercice  d’application: 

On  considère  deux  plaques  horizontales  parallèles,  de  grandes  surfaces,  séparées  par  une  couche  d’air, 
respectivement  aux  températures  T1  = 50°C  pour  la  plaque  supérieure  et  T2  = 30°C  pour  la  plaque 
inférieure.  La  distance  entre  les  deux  plaques  est  H = 5 mm.  L’émissivité  radiative  de  la  plaque 
chaude  est  prise  égale  à £1  = 0, 8 et  celle  de  la  plaque  froide  est  e2  = 0, 9. 


1.  Montrer  que,  dans  les  conditions  de  températures  considérées,  le  transfert  de  chaleur  par 
convection  naturelle  entre  l’air  et  les  plaques  ne  peut  pas  avoir  lieu.  Pour  cela,  calculer  le 
produit  du  nombre  de  Prandtl  par  le  nombre  de  Grashof  Pr  X Gr. 


> Pour  calculer  le  nombre  de  Nusselt  Nu,  on  utilise  la  corrélation  suivante  : 

Nu  = Cx(Prx  Gr)n 

On  rappelle  que  Gr  est  donné  par: 

p2  x g x p x AT  x H3 


Gr  = 


P 


ou  : 

A T=  \T1  - T2\ 





•hnH i a 
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P ( K x)  est  le  coefficient  de  dilatation  volumique  de  l’air  à pression  constante: 

1 

P = T~ 

1 m 

Tm  = — - — : température  moyenne  de  la  couche  d’air. 


Les  paramètres  C et  n prennent  respectivement  les  valeurs  : 

■ C = 0 et  n = 0 si  Pr  X Gr  < 3 X 105  -»  pas  de  convection ; 

■ C = 0,  27  et  n = si  3 X 105  < Pr  X Gr  < 3 X ÎO10  -*  le  régime  est  laminaire-, 

■ C = 0, 14  et  n = si  Pr  X Gr  > 3 X ÎO10  -♦  le  régime  est  turbulent. 


Fluide 

T t1  + t2 

2 

p (Kg/m3) 

Cp  (J/Kg.  °C) 

A (W/m.  °C) 

p ( Kg/m. s) 

Air 

40°  C 

1,136 

1004 

0,027 

0,0191x1 03 

2.  Calculer  la  densité  de  flux  de  chaleur  transféré  par  conduction  à travers  la  couche  d’air. 

3.  Calculer  la  densité  de  flux  de  chaleur  net  transféré  par  rayonnement  entre  les  deux  plaques. 

4.  En  déduire  la  densité  de  flux  de  chaleur  total  transféré. 

5.  Calculer  la  variation  relative  (en  %)  de  la  densité  de  flux  de  chaleur  total  si  l’on  utilise  un 
revêtement  sur  la  surface  froide  qui  ramène  son  émissivité  à e2  = 0, 1. 

On  donne:  cr=5,67x!0'8  W/m2.K4,  la  constante  de  Stefan-Boltzmann  et  g = 9, 81  m/s2, 
l’accélération  de  la  pesanteur. 


x Plaque  chaude  (S^) 

V Plaque  froide  (S2) 


Solution: 


1.  Montrons  qu’il  n’y  aura  pas  de  transfert  de  chaleur  par  convection  naturelle  entre  l’air  et  les 
deux  plaques. 


Pour  cela  calculons  le  nombre  de  Rayleigh  R a = Pr  X Gr: 
On  a : 

f 


p.  cp 

Pr  = ^r-~  0,71 


Gr 


p2  x g x p x AT  x H3 


277,2 


Ainsi  : 


Ra  = Pr  x Gr  = 196, 7 


Ra  = Pr  X Gr  = 196,7  < 3 X 105,  donc  il  n’y  aura  pas  de  transfert  de  chaleur  par  convection 
naturelle. 


2.  Densité  de  flux  de  chaleur  transféré  par  conduction  à travers  la  couche  d’air. 


<Pcd  = *x  — = 108  W/m 2 


H 


WV5  ^ 
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3.  Densité  de  flux  de  chaleur  net  transféré  par  rayonnement  entre  les  deux  plaques. 


avec: 


*Pr,net:l-*2  = ^12-  CT.  {T\  - T%) 


1 


Si 

+ ÇX 
s2 


Or:  F12  = 1 et  S±  = S2 
ainsi: 


Donc  : 


~rA «0,735 

1 1 

El  e2 

*Pr,net.'L—>2  ~ 102  W/m 


4.  Densité  de  flux  de  chaleur  total  transféré  : 


avec: 


tPtot  ^Pc v tPcd.  Vr.net, l->2 

A.  N.:  °cptot  -210  W/m2 

5.  Si  l’on  utilise  un  revêtement  sur  la  surface  froide  qui  ramène  son  émissivité  à e2  = 0, 1,  on  aura 

<P' r, net, 1^2  =Tl2 .0.(T\-T%) 

1 


Donc  : 


Pu—  T J ‘ 

1 — ; 1 

et  e2 


0,098 


<P'r, net, 1^2  * 13,  6 W/m 2 

La  densité  de  flux  de  chaleur  total  transféré  cp'tot  sera  égal  à: 

V tôt  Vcv  3r  Vcd  1"  V r.net, l->2 
^0 

A.  N.  : cp'tot  * 121,6  W/m2 

Il  en  résulte  que  la  variation,  en  pourcentage,  du  flux  de  chaleur  total  sera  égale  à: 

&<Ptot 

- -42% 

* P tôt 
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